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ABSTRACT 

The merger of two carbon-oxygen white dwarfs can lead to a spectacular transient - a type la supernova or 
accretion induced collapse to a neutron star - or the formation of a carbon star or massive, rapidly rotating white 
dwarf. Simulations of mergers have shown that the outcome strongly depends on whether the white dwarfs are 
similar or dissimilar in mass. In the similar-mass case, both white dwarfs merge fully and the remnant is hot 
throughout, while in the dissimilar case, the more massive, denser white dwarf remains cold and essentially 
intact, with the disrupted lower mass one wrapped around it in a hot envelope and disk. In order to determine 
what constitutes "similar in mass" and more generally how the properties of the merger remnant depend on the 
input masses, we simulated unsynchronized carbon-oxygen white dwarf mergers for a large range of masses 
using smoothed-particle hydrodynamics. We give full details of our results and tabulate key properties, in the 
hope of guiding analytical understanding and helping intercomparison between codes. Generally, we find that 
the properties of the remnants vary smoothly as a function of the two masses, with the remnant structure deter- 
mined primarily by the ratio of the central densities of the two white dwarfs. A density ratio of 0.6, equivalent 
to a mass difference of ^ 0.1 Mq, approximately separates similar and dissimilar mass mergers. Confirming 
previous work, we find that the temperatures of the merger remnants are not high enough to immediately ignite 
carbon fusion, except possibly for white dwarf masses approaching IMq. During subsequent viscous evo- 
lution, however, the interior will likely be compressed and heated as the disk accretes and the remnant spins 
down. From first-order estimates of the evolution, where we assume that the remnant spins down completely, 
that all rotational energy is used to expel matter to large distances, and that the remaining mass evolves adia- 
batically, we find that this can lead to ignition for many remnants. For similar-mass mergers, this would likely 
occur under sufficiently degenerate conditions that a thermonuclear runaway would ensue. 



1. INTRODUCTION 

A few percent of all white dwarfs (WDs) will eventually 
merge with another white dwarf. The outcome of such merg- 
ers will depend on the compositions of the WDs involved. For 
two helium WDs, a low-mass helium star might result, which 
would be observed as an sdOB star. For a helium WD merg- 
ing with a carbon-oxygen one, a helium giant could form, ob- 
servable as a hydrogen-deficient giant or R CrB star For two 
carbon-oxygen WDs (CO WDs), the outcome could vary be- 
tween simply a more massive WD, a carbon-burning star, an 
explosion, or collapse to a neutron star, depending on whether 
stable or unstable carbon fusion is ignited, and whether the 
total mass exceeds the critical mass for pycnonuclear ignition 
or electron captures (both close to the Chandrasekhar mass 
Mch). For mergers involving an oxygen-neon WD, the mass 
will always be high, and explosive demise or transmutation 
seems inevitable. 

The outcome of the merger of two CO WDs is uncertain in 
part because during the merger temperatures do not become 



hot enough to ignite sign ificant carbon fusion (e.g., [Loren- 
[Aguilar et al.| |2009| |LIG0 9 hereafter), except possibly for 
masses above ^0.9Mq (jPakmor et al. 2011, 2012). Hence, 
the final fate depends on subsequent evolutioriTin which dif- 
ferential rotation is dissipated, the remnant disk accretes, and 
the whole remnant possibly spins down. Due to these pro- 
cesses, the remnant could be compressed and heated, which, 
if it happens faster than the thermal timescale, would lead to 
increased temperatures and thus potentially to ignition. 

So far, efforts have focused on merging binaries with total 
mass M > Mch- The end result of such mergers is believed 
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to be either stable off-center carbon ignition, which would 
turn the merger remnant into an oxygen-neon WD a nd pos- 
sibly eventua lly result in accretion-induced collapse ( Saio &| 
|Nomoto|1998| ), or slow accretion, which allows the remnant to 
stay cool an d eventually ignite at high central density (Yoon] 
|et al.|[2007) l. Less massive mergers were usually thought to 
result in more massive, rapidly rotating CO WDs (Segretaui] 
et al. 1997; Kiilebi et al. 2010), but more recently it has been 
realized these might eventually become hot enough to ignite 



(van Kerkwijk et al. 2010 vKCJlO hereafter; IShen et al.|2012 
jSchwab et al.|2012| ). Indeed, vKCJlO argue that type la super 
novae result generally from mergers of CO WDs with similar 
masses, independent of whether or not their total mass ex- 
ceeds Mch (see below). For all these studies, the conclusions 
on whether and where ignition takes place depend critically 
on the structure of the merger remnnant. 

The merging process, and the merger remnant, have been 
studied quite extensively, mostly using smoothed-particle hy- 
drodynamics (SPH; e.g. Monaghan 2005 ). These simulations 
have shown that the outcome strongly depends on whether 
the WDs are similar or dissimilar in mass. In the similar-mass 
case, both WDs disrupt fully and the remnant is hot through- 
out, while in the dissimilar case, the more massive, denser 
WD remains essentially intact and relatively cold, with the 
disrupted lower mass one wrapped around it in a hot envelope 
and disk. Less clear, however, is what constitutes "similar- 
mass," and, more generally, how the merger remnant proper- 
ties depend on the initial conditions. 

In principle, for cold WDs of given composition, the rem- 
nant properties should depend mostly on the two WD masses, 
with a second-order effect due to rotation. In this paper, we try 
to determine these dependencies using high-resolution simu- 
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lations of WD mergers with the Gasoline SPH code, cover- 
ing the entire range of possible donor and accretor masses, 
but limiting ourselves to non-rotating WDs. Our primary 
aim is to identify trends between mergers of different masses, 
both to guide analytical understanding and to help scale other, 
perhaps more precise simulations. Here, our hope is that 
while the results of individual simulations may suffer from 
uncertainties related to the precise techniques and assump- 
tions used, the trends should be more robust. We also try 
to provide sufficient quantitative detail on the properties of 
merger remnants that it becomes possible to make analytical 
estimates or construct reasonable numerical approximations 
without having to run new simulations. 

Our work is co mplementary t o the rece nt surveys of rem - 
nant properties by 'Raskin et al. (2012 1 and Dan et al. ( 2012| l, 
in that they focus on different scientific questions (e.g., orbital 
stability; possible detonation). In contrast to our work, they 
assume that the WDs are co-rotating with the orbit. Whether 
this is a better assumption than no rotation depends on the 
strength of tidal dissipation, which unfo rtunately is not yet 
known (see Marsh et al. 2004; Fuller & Lai 2012). 

Our work also is part of a series of numerical studies inves- 
tigating the viability of sub-Chandrasekhar mass (sub-Mch) 
CO WD mergers producing SNe la, as proposed by vKCJlO 
This channel relies on similar-mass mergers producing rem- 
nants that are hottest near the center, and on compressional 
heating by subsequent accretion and/or magnetically medi- 
ated spin-down leading to ignition. The advantages of this 
channel are that it accounts for the absence of direct evidence 
for stellar companions, the observed SN la rate, and the de- 
pendence of SN la peak luminosity on the age of the host 
stellar population (because lower-mass merger constituents 
take longer to form). Since pure detonations of sub-Mch CO 
WDs produce li ght curves very similar to observed SNe la 
( [Shigeyama et al. 1992; Sim et al. 2010 ), it also removes the 
need for imposed deflagration-to-detonation transitions. Im- 
portant questions, however, remain, including what fraction 
of mergers leads to remnants that are hot near the center (in 
highly degenerate conditions), how the subsequent viscous 
phase proceeds in detail, whether ignition leads to a detona- 
tion, and whether the detonation of a remant that may still 
rotate and be surrounded by a disk would produce an event 
similar to an SN la. With our work, we attempt to address the 
first question. 

This paper is organised as follows. In Section [2] we de- 
scribe the SPH code we used, as well as our initial conditions. 
In Section|3] we present our results and give trends for a num- 
ber of pertment remnant properties. In Section |4] we test the 
robustness of our results, and in Section[5]compare our results 
with those of LIG09 and others. Lastly, m Section|6] we spec- 
ulate on the further evolution of ou r remnan ts, considering in 
particular whether, as suggested by |VKCJ10| some might lead 
to type la supernovae. 

2. CODE AND INPUT PHYSICS 

We simulate the mergers by placing non-rotating white 
dwarfs in a circular orbit with an initial separation aq chosen 
such that rapid mass transfer begins immediately. We then fol- 
low the merger for six orbits, at which time the remnant has 
become approximately axisymmetric. As in prior work, the 
morphology of all merger remnants is similar, consisting of a 
dense, primarily degeneracy-supported center surrounded by 
a partly th ermally-supported hot envelope (called a "corona" 
by |LIG09] ) and a thick, somewhat sub-Keplerian disk. We will 



use the terms "core", "envelope" and "disk" throughout this 
work. We also quite often refer to both the core and envelope 
simultaneously as the "core-envelope". 

2.1. The SPH Code 

With smoothed-particle hydrodynamics, one uses particles 
as a disordered set of interpolation points to determine contin- 
uum values of the fluid and model its dynamics. While gen- 
erally less accurate than grid-based codes, SPH often allows 
one to efficie ntly sim ulate phenomena with a large ra nge of 
lengthscales ( [Wadsley et al. 2004; P akmor et al.|2010| l. SPH 
is a Lagrangian method, meaning movement is automatically 
tracked, and regions of high density contain more particles 
and therefore are automatically more resolved. SPH has be- 
come the method of choice for merger simulations, and partly 
because we wish to compare with previous work, we chose 
SPH as weU. 

For our simulations, we use Gasoline ( [Wadsley et al.|20()4l ), 
a modular tree-based SPH code that was designed and has 
been used for a wide range of astrophysical scenarios, from 
galaxy interactions to planet formation. It aims for tight con- 
trols on force accuracy and integration errors. In our sim- 
ulations, total energy is on average conserved to 0.3%, and 
angular momentum to 0.006%. 

By default. Gasoline uses the standard Monaghan and Gin - 
gold formulation for artificial viscosity (see |Monaghan|2005 1, 
augmented with a Balsara switch to reduce viscosity in non- 
shocking, shearing flows. IGuerrero et al.|([2004) found that 
such a prescription did not reduce viscosity sufficiently, re- 
sulting in exc ess spin-up of the remnant core and associated 
shear heating. |Yoon et ar| ( |2007| l used variable coefficients for 
the linear and quadratic viscosity terms in the SPH equations 
of motion and energy, setting these values to a = 0.05 and 
/3 =0.1, respectively, where shocks are absent, and around 
un ity where they are present. A similar formulation was used 
in |Dan et al.| ( |20lT] l. Since Gasoline includes it as well, we 
have used it for our study. Excess viscosity nevertheless re- 
mainsa potential problem; we investigate its effects further in 
Sec. 1431 

We modified Gasoline to include support for degenerate gas 
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through the H eknholtz equation of state (EOSV ([Timmes & 

Irnno] 



|Swesty|2000[ ). This code interpolates the Helniholtz free en 
ergy of the electron-positron plasma, along with analytical ex- 
pressions for ions and photons, to determine pressure, energy 
and other properties from density and temperature. It is fast, 
spans a large range of density and temperature, and has, by 
construction, perfect thermodynamic consistency. To obtain 
quantities as a function of density and internal energy, as used 
in Gasoline, we utilized a Newton-Raphson inverter, disabling 
Coulomb corrections to keep the energy-temperature relation 
positive-definite. 

A problem we encountered in modelling degenerate mate- 
rial is that it is affected much more strongly by a ubiquitous 
noise source in SPH simulations: particles moving ran domly 
at a small fraction of the speed of sound (e.g., Springel 2010[ ). 
Artificial viscosity helps damp this noise, but generates cor- 
responding thermal energy. This spurious thermal energy is 
typically a small fraction of the internal energy, and therefore 
if the pressure support is predominantly from thermal pres- 
sure changes to the temperature are negligible. In our sys- 
tems, however, degeneracy pressure dominates, and the small 
changes in internal energy correspond to large temperature 
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changes, at times comparable to the physically expected val- 
ues. We characterize this "spurious heating" in Sec. 4.6 and 
show that it does not unduly affect our work's conclusions. 
However, it makes it difficult to run much longer simulations. 

Another issue we encountered is that a particle's energy can 
drift below the Fermi energy. To prevent pressure from be- 
coming unphysical, we do not let it decrease below its zero- 
temperature value. We leave the energy as is to maintain en- 
ergy conservation, except we do not let it decrease below half 
the Fermi energy. The latter is to prevent particle energies 
from reaching below zero, which under rare circumstances 
appears to happen when particles perform a great deal of me- 
chanical work. We found this occurs primarily for particles 
that are flung out of the system by the merger and are cooling 
rapidly, and therefore are confident it has only a very minor 
effect on our simulations. 

In our work, we ignore outer hydrogen and helium layers, 
composition gradients, and any nuclear reactions. This is pri- 
marily because previous work has found that nuclear process- 



ing was unimportant during the merger. For instance, |LIG09 
found fusion released ^ 10'*' erg for their 0.6 - 0.8 merger, 
orders of magnitude smaller than the ^ 10^" erg binding en- 
ergy of the remnant. Only for mergers involving very mas- 
sive, >O.9M0 WDs might this assumption bre ak down, with 
the possibility of carbon detonations arising (Pakmor et al.' 
ITOlO 2011 201 2) but see|Raskin et al. 2012; Dan et al. 2012). 
Similarly, Raskin et al. (2012 1, who included standard helium 
envelopes of ^ 1 -2% of the WD mass in their simulations, 
found that only for accretors with masses above ^ 1 Mq did it 
made a substantial difference: a helium detonation would in- 
ject ^ lO'*^ erg into the merger remnant. While this led to addi- 
tional heating, it was insufficient to trigger substantial carbon 
burning or unbind any portion of the remnant (helium deto- 
nations have also bee n found for lowe r-mass accretors with 
co-He hybrid donors ; |Dan et al |20 1 2) . 

2.2. Input Stars 

We created spherical white dwarfs using pre-relaxed cells 
of particles rescaled to follow the appropriate enclosed mass- 
radius relation determined using the Helmholtz equation of 
state. We assumed a composition of 50% carbon and 50% 
oxygen by mass, and a uniform temperature of 5 x lO^K. 
The stars were then relaxed in Gasoline for 81 s with ther- 
mal energy and motion damped (to 5 x 10''K and cm s~', 
respectively) during the first 41s. Ultimately, however, par- 
ticle velocity noise (described above) prevented cooling of 
> 5 X 10^ gem"-' material to below 10^ K. We checked that 
after relaxation the density profile was consistent with the so- 
lution from hydrostatic equilibrium. We found that the central 
densities agreed to within 2%, but the radii, as defined by the 
outermost particle of a relaxed WD, were on average about 
7% too small, reflecting our inability to model the tenuous 
WD outer layer^ 

We used a constant particle mass of 10^** g, so that a O.4M0 
WD has 8 X lO"* particles, and a I.OM0 WD has 2 x 10^ 



by |Dan et al. (2012). [Raskin et aLl d20T2| performed a 
olution test for a merger of two 0.8 IM^ WDs, varying 



These numb ers are similar to those used by LIG09 and Yoon 
et al. (20071, and exceed the ~ 2x10* part icles per star used 

res- 
ying the 

^ Our relaxed WD.s also show evidence of sub-kernel radial banding of 
particles, which does not appear in any interpolated quantities. We do not 
believe this banding has an effect on our simulations except for a possible 
reduction in effective resolution, but will investigate remedies in future work. 



number of particles per star from 10^ to 2 x 10*. They found 
differences of ^ 2% in the mass of the core plus envelope, 
disk half-mass radius, and inner disk rotation frequency. The 
one qualitative difference they found was that at their highest 
particle resolution, the WDs failed to break symmetry and dis- 
rupt (note that they assumed co-rotating WDs, making such a 
stable cont act c onfiguration possible). We perform our own 
test in Sec. 4.7 and find similar results, making us confident 
that our resolution is sufficient. 

2.3. Input Binaries 

We relaxed 0.4, 0.5, 0.55, 0.6, 0.65, 0.7, 0.8, 0.9 and 1.0 Mq 
white dwarfs, and combined them in all possible permutations 
to form our parameter space of binaries. These values were 
chosen to represent the range of possible CO WD masses, 
with greater resolution near the empirical peak at ^0.65Mq 
of the mass distribution of (single) CO WDs (Tremblay &j 
Bergeron 2009) . We also performed additional simulations 
with 0.5'/5 - 0.'65, 0.625 - 0.65 and 0.64 - 0.65 M© binaries 
to explore the outcomes of similar-mass mergers. We thus 
simulated 48 mergers in total. 

We placed two relaxed, irrotational WDs in a circular orbit. 
We chose the initial separation «(> such that the donor WD just 
fills its Roche lobe, taking the location of the donor's outer- 
most particle as its radius and using the Roche lobe approxi- 
mation (for a synchronized binary) from Eggleton ( 198 3|l. 
This simple initial condition is similar to that of Pernor 



et al. (2010), and implies that the binary system as a whole is 
not equilibrated. Therefore, as the simulation begins, the two 
WDs react to the tides, become stretched, and strong Roche 
lobe overflow ensues because the donor overshoots its Roche 
radius (in a widely separated binary, the donor would start to 
pulsate). As a result, the donor disr upts after j ust on e to two 
orbits. For synchronized binaries, Dan et aL] ( |20i r) showed 
that the onset of mass transfer is much more gentle if the WDs 
are relaxed in the binary potential, disruption occurring only 
after several dozen orbital periods. They also showed that 
this results in systematic changes in the merger remnants. It 
is not clear whether the same will hold for unsynchronized 
binaries, since the accretion stream hits a surface that, in its 
frame, counterrotates, and therefore accretion is always much 
less gentle than for synchronized WDs. The difference is par- 
ticularly dramatic for similar-mass binaries, where, in the syn- 
chronized case, the WDs can come into gentle contact, while 
in the unsynchronized case, any contact is violent. Unfortu- 
nately, it is difficult to test the effect of proper equilibration for 
unsynchronized binaries, since one has to relax to non-trivial 
initial c ond itions. A better approxi mation was attempted by 
|LIG09| and [Guerrero eCal] ^M^, who started their WDs 
further out and reduced the separation artificially until mass 
transfer began. In their simulations, disruption still followed 
very quickly. Given that, and wanting to avoid any partial 
synchronization, we kept our simpler setup, and tested it by 
running sim ulati ons with varying ao. We will discuss these 
tests in Sec. 4.2 and compare our results with those of others 
inSeclSl 

2.4. Merger Completion Time 

It is difficult to decide when a merger is "complete", since 
for some cases remnant properties continue to evolve long 
after the two WDs coalesce, with (artificial) viscosity redis- 
tributing angular momentum and heating the remnant. As a 
visually inspired criterion, we decided initially to use the de- 
gree of non-axisymmetry, continuing simulations until they 
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were less than 2.5% non-axisymmetric, as measured from the 
ratio of zeroth to largest non-zero Fourier coefficient of par- 
ticles binned in azimuth. However, this had its own issues: 
in dissimilar-mass mergers - where most of the particles are 
in the accretor, already roughly axisymmetric following the 
merger - our convergence criterion was achieved while the 
outer disk was still obviously non-axisymmetric. In equal- 
mass mergers, which are inherently more axisymmetric, com- 
pletion also was too soon, before the densest material had 
reached the center of the remnant. 

For the majority of our systems, however, the time required 
to reach 2.5% non-axisymmetry was roughly constant in units 
of the initial orbital period, at 6.1 ± 1.2. For about the same 
time, axisymmetry was also achieved (by subjective visual 
inspection) for both dissimilar-mass mergers (except, in ex- 
treme dissimilar-mass cases, the outermost regions of their 
disks) and for equal-mass mergers (where the densest material 
had reached the center). We therefore use 6 orbital periods of 
the initial binary as the completion time of our simulations. In 
Sec. 4.4 we discuss the effect of continuing our simulations 
for 2 further orbital periods. 

3. RESULTS 

With our 48 simulated mergers in hand, we try to determine 
scaling relations of global quantities such as the remnant and 
disk mass, highest temperature, etc., and look for homologies 
in the remnant profiles. For our analysis, we use a cylindri- 
cal (C3,0,z) coordinate system centered on the remnant core. 
Properties on the equatorial ((53, 0) plane - defined as the orig- 
inal orbital plane - are averaged over </> using particles within 
jh^ of the equatorial plane, where /jz is the remnant's rota- 
tional axis (to = 0) central scaleheight (see Sec. |3.2.2"l l. Proper- 
ties along the rotational (z) axis are averaged within a cylinder 
W < j/iz. We use -j^/iz as the bin size along both the equatorial 
plane and rotational axis. We determine properties mostly as a 
function of enclosed mass M(r), which we define sphericall}]^ 
Thus, we show, e.g., equatorial plane temperature r(!I!) as a 
function of M(r = 03), the mass enclosed within a sphere with 
radius r = 03. 

3.1. Representative Mergers 

As found for previous simulations, qualitatively the most 
important factor controlling the merger outcome is whether 
the WD masses are "dissimilar" or " similar". In the for- 
mer case, where the donor is significantly less massive than 
the accretor, only the donor overflows its Roche lobeljis dis- 
rupted, and accretes onto the accretor. The accretea mate- 
rial is shock-heated on impact, lifting degeneracy. Hence, the 
merger remnant consists of a non-degenerate hot envelope and 
small, thick Keplerian disk, both surrounding a cold core con- 
taining the largely unaffected accretor. 

In the latter case of a similar-mass merger, there is a large 
degree of mixing between the two stars. For exactly equal 
masses, both stars are disrupted simultaneously, and their ac- 
cretion streams impact each other near the system's barycen- 
ter Material from the centers of both stars initially forms 
a thick, cold, dense torus orbiting the barycenter; this torus 
slowly shrinks due to viscous drag, pushing the shock-heated 

^ Arguably, enclosed mass is more properly defined within equipoten- 
tial surfaces, but this makes comparison with other simulations harder. For 
dissimilar-mass mergers, the difference is slight. 

The lower mass WD is larger and thus always fills its (smaller) Roche 
lobe first. 



accretion stream material above and below the equatorial 
plane. When the stars have slightly different mass, the lower- 
mass one disrupts first, forming an accretion stream (or series 
of streams) that penetrates into the center of the accretor (re- 
gardless of whether or not the other also disrupts). 

We show the differences between similar and dissimilar- 
mass merger remnants using two representative examples in 
Figs. [T]and|^ a 0.4 - 0.8 Mq highly dissimilar and a 0.6 - 
0.6 Mq equal-mass merger, respectively. One sees that the 
remnant morphologies are very different, consistent with pre- 
vious work. The 0.4 - 0.8 Mq merger features a cold, nearly 
non-rotating and thus spherically symmetric remnant core, 
surrounded by a hot envelope with roughly equal degeneracy 
and thermal support, which itself is surrounded on the equa- 
torial plane by a rotationally supported non-degenerate thick 
disk that holds most of the angular momentum. The accretor 
forms the core, largely undisturbed by the merger, while the 
envelope and disk are composed almost entirely out of donor 
material. The hottest points are on the interface between the 
core and the envelope]^ The 0.6 - 0.6 Mq remnant, on the 
other hand, has a massive, warm, partly rotationally supported 
and thus ellipsoidal core, and a very small but thick disk, both 
of which consist of material from both stars. No distinct en- 
velope is formed. The hottest points are within the remnant 
core, just above and below the equatorial plane, arising from 
accretion stream material pushed out by the shrinking dense 
torus. 

A good way to visualize how mergers transition between 
dissimilar and similar-mass is to look at changes in the rem- 
nant properties with varying donor mass. In Fig. [3] we show 
curves for accretors of 0.65 (left) and l.OM0(right). One sees 
that remnants of highly dissimilar-mass mergers, with mass 
ratio = M^/M-j, < 0.5, have properties resembling the 0.4 - 
0.8 Mq merger: their donor and accretor barely mixed, their 
temperature curves have off-center hot plateaus, and their an- 
gular velocity profiles feature an off-center bump. The equal- 
mass, = 1 cases resemble the 0.6 - 0.6 Mq remnant: they 
have flat temperature profiles and centrally peaked angular ve- 
locity profiles. Intermediate cases have intermediate profiles, 
with the bumps in the temperature and angular velocity pro- 
files widening with increasing q^. The 0.4 - 0.8 Mq and 0.6 - 
0.6 Mq remnants therefore lie at the extremes of what merger 
remnants look like. 

The similarity between some of the curves for the 0.65 and 
1.0 M0 accretors in Fig. [3] suggests a homology. The sim- 
ilarity is closest for mergers with the same mass difference 
AM, as can be seen in Fig. [4] For equal-mass mergers, all 
profiles are similar, simply scaled by a factor that depends 
on the total mass (except the 1.0 - 1.0 Mq merger; see be- 
low). As AM increases, the profiles are slightly less simi- 
lar: with increasing total binary mass, the degree of mixing 
decreases, and the temperature and angular velocity maxima 
drift to slightly lower fractional enclosed mass. Nevertheless, 
the profiles still resemble one another far more closely than 
they resemble curves with other AM. The same holds for 
profiles along the rotational axis. 

It may seem surprising that the controlling parameter be- 
tween these approximate homologies is the mass difference 
AM rather than the mass ratio q^. Empirically, however, the 
case is clear: e.g., the 0.4 - 0.5 (second column, yellow) and 
0.8 - 1.0 Mq (third column, black) mergers have the same 
qm, but different AM, and their structures clearly differ from 



' The higher temperatures near the core are spurious; see Sec. 4.6 
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Figure 1. Structure of a 0.4 - 0.8 Mq merger remnant, representing the general outcome of a merger of white dwarfs with dissimilar mass. Upper left and middle 
- maps of density p and temperature T along slices in the xy and jcz-planes. Lower left - maps and contours of density, temperature, and angular frequency Q, in 
the (05, z) plane, averaged over cylindrical coordinate and over itz (with 1 added to f2 to avoid problems with the logarithmic intensity scale). Middle - enclosed 
masses of donor and accretor material and Ma (solid red and blue, resp.), and fraction of donor material /d at a particular mass shell (dashed magenta). Middle, 
one but lowest - temperature-density profile with enclosed masses in 0.2 Mq increments indicated, both along the equatorial plane (solid curve, squares) and 
along the rotational axis (dot-dashed curve, circles). Middle, bottom - enclosed mass as a function of r, with the total mass indicated by the horizontal dashed red 
line. Right-hand column, top to bottom - density, temperature, entropy, angular (cyan) and Keplerian (blue) frequency, and degeneracy (blue), thermal (red) and 
rotational (cyan) energy densities as a function of enclosed mass M, both along the equatorial plane and along the rotational axis (solid and dot-dashed curves, 
respectively). In all graphs, the start of the disk (where the centrifugal acceleration equals half the gravitational one) and the mass enclosed within the equatorial 
radius with maximum temperature are marked by vertical green and blue dashed lines, respectively. [See the electronic edition of the Journal for Figs. 1.1-1.48.] 
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Figure 2. As Fig. Set[T] but for a 0.6 - 0.6 Mq merger remnant, representing tlie general outcome of a similar-mass merger. 
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one another. The same is true for the 0.4 - 0.6 (third column, 
cyan) and 0.6 - 0.9 (fourth column, brown) Mq mergers. As 
we discuss below, the similarity of mergers of similar AM 
likely reflects the close relation between the ratio of central 
densities and mass difference. 

Before discussing the homologies and trends further, we 
should note the one dramatic exception. The 1.0 - 1.0 Mq 
simulation differs fundamentally from its fellow AM = 
mergers. During the evolution of this system, unlike for all 
other equal-mass mergers, one WD was disrupted long before 
the other, and as a result material from one star (arbitrarily 
designated the "donor" before the start of simulation, hence 
the "inverted" mixing profile in Figs. [3] and |4]) preferentially 
resides near the center of the remnant. This system also of- 
ten appears as an outlier in Sec. |3.2| below. The 0.9 - 0.9 Mq 
merger also did not have equal mixmg b etween the two stars, 
though the difference is much smaller. |Raskin et al.| ( |2012] l 
noticed the same effect in their simulations, and concluded it 
reflected the fact that more massive WDs are much more con- 
centrated and therefore harder to disrupt. This seems a likely 
explanation. 

3.2. Merger Trends 

A major goal of our work is to establish how various 
global properties of the merger remnant, such as remnant 
core and disk mass, maximum temperature, and maximum 
angular velocity, vary as a function of accretor and donor 
mass. By quantifying these trends, we hope to help develop 
a parametrized model of merger remnants. Before discussing 
trends, however, we stress that they are necessarily approxi- 
mate - second order effects, numerical noise and our choice 
of stopping time all affect the remnant properties. Moreover, 
while integrated values like total thermal energy do not fluctu- 
ate from timestep to timestep, values at specific points in the 
remnant do (as noted the following sections). For instance, the 
mass enclosed within the radius of peak equatorial temper- 
ature becomes ill-defined for similar-mass mergers because 
these have rather flat temperature profiles (Fig. [3]|. To partly 
mitigate these fluctuations, the values presented below were 
determined by averaging frames from the simulation over an 
eight second span, centered on the time corresponding to six 
orbits of the initial binary. 

As might be expected from the approximate homologies 
described above, we found that many properties scaled well 
with AM. Of course, a scaling with a dimensional mass dif- 
ference makes little sense; we believe its success reflects the 
fact that over the range of 0.4 - I.OMq, the central density 
Pc depends approximately exponentially on mass, with pc ~ 
3.3 X lO^gcm-^exp[5.64(M/M0-l)] (see Fig.^. Hence, a 
given mass difference AM corresponds to a given ratio of cen- 

Pc,d/Pc,a has 



0.65 M„ Accretor 1.0 M,^ Accretor 



3.2.3 



tral densities, Pcd/Pca- As argued in Sec. 
a straightforward interpretation: it characterizes the degree of 
mixing between the donor and accretor We therefore discuss 
trends as a function of qp = /9c,d/Pc,a from hereon. Where 
necessary, we refer to the mass ratio as q^^. 

3.2.1. What Constitutes Similar-Mass? 

As qp increases from a small value toward unity, the 
merger remnant's morphology shifts from resembling Fig. [T] 
(dissimilar-mass) to resembling Fig. [2] (equal-mass). From 
Fig. [3j ones sees that there is no particular qp at which one 
transitions from "dissimilar" to "similar." Nevertheless, we 
can determine a rough critical value of qp that separates merg- 
ers in which the core is largely unaffected from those in which 




0.4 0.6 



0.4 0.6 



Figure 3. Properties of mergers with 0.65 M0 (left) and I.OMq (right) ac- 
cretors, for donor masses of 0.4 (red), 0.5 (orange), 0.55 (lime), 0.575 (green), 
0.6 (cyan), 0.625 (light blue), 0.64 (blue), 0.65 (dark blue), 0.7 (magenta), 0.8 
(purple), 0.9 (brown), and I.OMq (black). Shown are, from top to bottom, 
density p, fraction of donor material /j, angular frequency Q, temperature T, 
thermal energy iJth. ™d rotational energy £rot, all as a function of fractional 
enclosed mass M/Mtot- All properties are determined along the equatorial 
plane, except for /j which is defined spherically. The 1.0- 1.0 Mq merger 
(dashed black line) is an outlier; see text. 

it is changed significantly, a separation that likely affects the 
outcome of post-merger evolution. 

To determine the critical value, we show in Fig.[6]the ratio 
of central to maximum temperature, Ti/rmax, central to maxi- 
mum angular velocity, ftc/^max, and the fraction of donor to 
accretor material within the central core, (/d//a)cc, where we 
define the central core as a sphere with radius /iz. All three 
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AM= 

Pc.d/Pc.a = 1 



AM = 0.1 

Pcd/Pca = 0.57 



AM = 0.2 

Pcd/Pca = 0.32 



AAf = 0.3 
Pcd/Pca =0.19 




0.0 0.2 0.4 0.6 

M/M,,, 



0.0 0.2 0.4 0.6 

M/M,„, 



0.0 0.2 0.4 0.6 0.1 



0.0 0.2 0.4 0.6 0.8 1.0 

M/M,,, 



Figure 4. Dependence of the properties of mergers on mass difference, with, from left to right, AM = Ma 
shown, coloring, and line styles are as in Fig. [3] except color represents accretor mass. 



M(i = 0.0, 0.1, 0.2, and 0.3 Mq mergers. Properties 



properties are measures of the extent to which the core has 
been affected: mixed regions tend to be hotter and more spun 
up, and contain material from both stars. 

From Fig. [6] one sees that flc/^mux approaches unity at 
qp ~ 0.6; at nigher values, the angular velocity profile has 
a plateau or central peak rather than an off-center bump. Also 
at qp ~ 0.6, (/d//a)cc starts to deviate from zero, i.e., donor 
material begins to penetrate the central core. The tempera- 
ture points show the transition is not abrupt: T^^/T^ax starts 
to deviate from its downward trend (which reflects spuri- 
ous heating in the most dissimilar-mass mergers; Sec 
around qp ~ 0.3 and continues to increase until qp 



4^1 
l.OTat 



qp ~ 0.6, Tc/Tinax — 0.5. Overall, this suggests that while the 
dependence is gradual, the morphology changes most around 
qp ~ 0.6. This conclusion is confirmed by looking at the two- 
dimensional remnant temperature structures (Figs. [I] and |2]). 
At qp <C 0.6, the remnant core has a large, spherically sym- 
metric cold region, the nearly unperturbed accretor. This cold 
region shrinks with increasing qp, and at qp ~ 0.6, spherical 
symmetry is broken. For still larger qp, the cold region be- 
comes a flat slice sandwiched between hotspots off the equa- 
torial plane. 

Given the above, we define "similar-mass" mergers as those 
with donor to accretor central density ratio qp > 0.6, and 
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0.6 0.1 

M (M„) 

Figure 5. Relation between central density pc and mass M for white dwarfs, 
showing both the results of relaxing white dwarf models in Gasoline (red 
points), and integrating hydrostatic equilibrium directly for spherically sym- 
metric, non-rotating carbon-oxygen white dwaii's with T = 5 X 10*K (blue 
line). For the mass range considered, the central density depends roughly 
exponentially on mass. 
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Figure 6. Dependence of merger core properties on the ratio of the donor 
and accretor central densities, Pc.d/Pc.a. Shown are the ratio of central to 
maximum temperature Tc/7max (squares), central to maximum angular ve- 
locity sic/n 

max (triangles), and central core donor to accretor mass frac- 
tion (/d//a)cc (circles), with colors representing different accretor masses, 
encoded as in Fig. |3] The vertical line marks qp = Pc.d/Pca = 0.6, where 
f2c/f^raax reaches unity, (/d//a)cc becomes non-zero, and Tc/Tm^j. ~ 0.5. We 
suggest it separates "dissimilar" from "similar" mass mergers. 

"dissimilar-mass" mergers as those with qp < 0.6. This crit- 
ical density ratio corresponds to a mass difference AM ~ 
O.IMq. 

3.2.2. Structural Trends 

We begin our discussion of trends with size and density pa- 
rameters. 

The rotational axis central scaleheight. — We define the rota- 
tional axis central scaleheight as the characteristic width a 
of a Gaussian fit to the density distribution along the z axis at 



(53 = 0. /iz is a measure of the vertical extent of the remnant. 
We find that the ratio h^/h^, where is the central scaleheight 
of the accretor, is reasonably well-approximated by. 



1.03 



(±0.02), 



(1) 



where the uncertainty Usted in parentheses represents the root- 
mean-square (RMS) of the residuals around the approxima- 
tion (see Fig. [7] a).). For highly dissimilar-mass mergers, 
approximately equals the scaleheight of the accretor, while for 
similar-mass mergers, /iz is lower due to rotational support. 

The vertical scaleheight increases with increasing (D: the 
scaleheight at the location of maximum temperature, /i(i;nax), 
ranges from to 1.21/iz, and the scaleheight at maximum 
angular velocity /i(ilmax) ranges from to 1.88/;z. The pref- 
actor for both heights increases with increasing accretor mass 
Ma and decreasing qp. 

The equatorial plane central scaleheight. — Similar to /iz we de- 
fine /in - the characteristic width of a Gaussian fit to the den- 
sity distribution along the equatorial plane - as a measure of 
the equatorial extent of the remnant. The ratio hm/ha can be 
parametrized by 



hm 



= 0.96-1-0.89^2 (±0.08), 



(2) 



where we excluded the 1.0 - 1.0 Mq merger remnant for our 
fit (see Fig. l7]a).). The dependence on increasing qp reflects 
the increaseorotational support of the remnant core. 

The central density of the remnant. — The central density, p^, is 
always within a factor two of the central density of the accre- 
tor, pc.a- In Fig. [7|b)., one sees that for given accretor mass, 
pjpc.a increases with increasing qp for highly dissimilar-mass 
mergers due to increasing compression of the remnant core, 
but begins to decrease because of increasing rotational sup- 
port around qp ~ 0.3. We could not find a simple parametriza- 
tion for these curves. Note that for some systems, pjpc.n con- 
tinues t o inc rease as we continue our simulation. As discussed 
in Sec. |2.4| this is probably because artificial viscosity forces 
the merger remnant to undergo accelerated viscous evolution. 

3.2.3. Mass Distributions 

The merger mixes material between the donor and accretor. 
Here, we describe how this changes as a function of qp, as 
well as how the material is distributed between the pressure- 
supported core and envelope and the rotationally supported 
disk. 

The masses of the core-envelope and disk. — We formally de- 
fine the core-envelope as the part of the remnant inside the 
inner disk radius, i.e., that is supported primarily by pressure 
(degeneracy for the core, thermal for the envelope) and not 
rotation. Since in every merger very little mass is ejected, a 
trend for either the core-envelope or the disk mass (Mce and 
Mfjisi;, resp.) suffices to determine both. The ratio of M^e to 
the accretor mass Ma is well described by. 



Mce 

M„ 



l-l-0.81^p 



(±0.03), 



(3) 



if the 1.0 - 1.0 M0 merger is neglected, and the fit's y-intercept 
is forced to unity. See Fig.|8]a). 
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Figure 7. Structural properties of mergers, (a) Central scaleheights along 
the rotational axis (circles) and along the equatorial plane (squares), scaled 
to the scaleheight of the accretor, ftz/^a and /im/'ia- (b) Central density of the 
merger remnant scaled to the central density of the accretor, Pc/Pc,a- Colors 
represent different accretor masses, encoded as in Fig.|3] Triangles represent 
the outlying 1.0- 1.0 Mq merger. 



The mass enclosing 50% of the donor material. — The further 
the donor penetrates, the smaller will be the mass enclosing 
half the donor's material, MgncCjMd). For mergers with qp < 
0.8, Menc(5Md)/Ma ~ 1.30, with an RMS residual of 0.03 
(Fig.[8]c).). We present this trend mostly because we discuss 
similar thermodynamic and rotational enclosed masses, but it 
is somewhat difficult to interpret physically, since Menc( jA^d) 
increases with donor mass but decreases with mixing, which 
also depends on donor mass. The trend is easier to interpret 
using enclosed accretor mass rather enclosed total mass, as 
done below. 

The accretor mass enclosing 50% of the donor material. — As a 
different measure of the depth to which the donor penetrates, 
we consider just the accretor material within the mass enclos- 
ing half the donor, Mgnd^Md)- jMd. This should equal the 
accretor mass if the donor is deposited above the accretor, and 
half the accretor mass if the two stars are completely mixed. 
For qp < 0.8, it can be approximated by (Fig.lslb).), 



Me„c(iMd)- 



iMd 



where we forced the intercept to be unity. In this regime, 
roughly half of the donor remains outside of the accretor, 
though the trend discussed next indicates that the other half 
which does penetrate the accretor is spread across a much 
larger region at higher qp. When qp > 0.8, the ratio drops 
sharply downward, indicative of the more thorough mixing 
expected for the similar-mass case. However, the existence 
and exact location o f thi s drop may be a function of initial 
conditions (see Sect. |4.2| l. 

The region over which the donor is spread. — As a measure of the 
thickness of the region affected by the merger, we use the dif- 
ference of the mass enclosing 75% of the donor material with 
that enclosing 25% of the donor material, i.e., AMenc(A^d) = 
^enc( jA/d)-Menc( jMd). Since 50% of the donor is within this 
range, AMenc(A^d)~ |-^d is a measure of the amount of accre- 
tor mixed with the donor. For qp < 0.8, the ratio of the latter 
to the total accretor mass follows. 



AMe„c(Md)- 



iMd 



= 0.30qp (±0.02), 



(5) 



while for qp > 0.8 the trend curls upward until it reaches 0.5, 
the value expected for completely mixed remnants. See Fig. 
id). 

Combining the two above trends, we can formulate a qual- 
itative picture of mixing. For qp < 0.8, the donor can be 
thought of as being deposited onto the accretor and mixing 
with the accretor's outer layers, while for qp > 0.8, the accre- 
tor also disrupts substantially, leading to a regime where both 
stars mix more uniformly. The region over which the donor 
is spread, or thickness of the mixed layer, in both cases de- 
pends on qp, which suggests that the relative densities of the 
donor and accretor govern mixing, i.e., the donor mixes sig- 
nificantly with all accretor material up to some fraction of the 
central density of the donor. Additional evidence of this will 
be seen in the thermodynamic trends below. 

One might consider an alternate picture in which the donor 
dredges up a constant fraction of its own mass in accretor ma- 
terial. If this were the case, we would expect (AMenc(Md)- 
jMi)/Md to roughly be constant. Our results, however, show 
that this quantity is nearly a straight line that is close to zero 
for small qp and ^ 0.5 when ~ 1. This seems more consis- 
tent with mixing being determined by density. 

3.2.4. Energy Balance 

The energy balance of the remnants indicate their primary 
means of support. Since the remnants are virialized, we con- 
sider how the ratio of degeneracy, thermal, and rotational en- 
ergy to the total internal energy of the remnants varies with 

Energy balance of the entire remnant. — The support against 
gravity changes from being due mostly to degeneracy pres- 
sure at low qp to having a substantial rotational contribution 
at q'p ~ 1 . This is because for highly dissimilar-mass mergers 
most of the internal energy is locked up within the accretor, 
which is hardly heated or spun up. For similar-mass mergers, 
however, the donor penetrates into the center of the accretor, 
causing heating and spin-up. 

The total gravitational potential energy of the merger rem- 
nant can be described adequately by a constant fraction of 



= l-0.190^p (±0.009), (4) 



■^pot 



GMIJR.^ 



= 0.49 (±0.01). 



(6) 
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Figure 8. Mixing, iieating, and spin-up (left to right) for mergers, (a) Scaled mass of the remnant core-envelope (where scaling here and below is to the 
accretor mass), (b) Fraction of the accretor within the mass enclosing half the donor mass, (c) Scaled remnant mass enclosing half the donor mass, (d) 
Fraction of the accretor mass within the region enclosing 25-75% of the donor mass, (e) Maximum equatorial temperature Tmax (circles), with the approximation 
Tinax = 0.20GMiimp//:BRa overdrawn. Maximum temperatures along the z-axis are shown with crosses, (f) Density at the location of Tmax (symbols as above), 
(g) Scaled mass enclosed within the radius of Tmax (symbols as above), (h) Scaled mass enclosing half of the remnant thermal energy, (i) Scaled mass of the 
region enclosing 25 - 75% of the remnant thermal energy, (j) Maximum angular velocity fimax (circles) with best fit Qmax = S.SHorb overdrawn, (k) Fraction 
of the angular momentum in the core-envelope. (1) Scaled mass enclosed within the radius of maximum angular velocity, (m) Scaled mass enclosing half of the 
total remnant rotational energy, (n) Scaled mass of the region enclosing 25 - 75% of the remnant rotational energy. Colors represent different accretor masses, 
encoded as in Fig. [5] Triangles represent equatorial plane values, and x-marks rotational axis values, of the 1.0 - 1.0 Mq merger. 
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From the virial theorem, the internal energy should be re- 
lated to the potential energy by 3((7) -l)Ei = --Epot, where (7) 
is an appropriately averaged equivalent to the adiabatic index. 
Since our remnants have cores where the electrons are be- 
coming relativistic, one has (7) somewhat smaller than 5/3, 
especially for the more massive remnants. We find that the 
ratio Ei/\Epot\ can be described by. 



"pot I 



M 

:0.18 — + 0.42 



(±0.01), 



(7) 



which is between ^0.5 and 0.6 for low and high Ma, respec- 
tively. 

The fraction of the internal energy carried by degeneracy 
and rotation is fairly well described by. 



^10 



= 0.31^ 



1/2 



(±0.01), 



= 0.92-0.34^ 



1/2 
p 



(±0.02). 



(8) 
(9) 



With these, the fraction carried by thermal energy can also 
be calculated; as shown in Fig. |9] a)., the fraction in ther- 
mal energy first increases with increasing qp, but turns over 
at qp ~ 0.7, decreasing afterwards. This reflects the compe- 
tition between increased thermal energy from the two stars 
mixing, and increased rotational support from the spin-up of 
the core. 

Overall, for highly dissimilar-mass mergers, the internal 
energy is partitioned into degeneracy, rotational and thermal 
energy with a ratio of approximately 8:1:1, reflecting that, 
as stated above, such mergers are almost entirely supported 
by degeneracy pressure. Similar-mass mergers, on the other 
hand, partition their internal energies with the ratio 6:3: 1, i.e., 
rotational support is significant. 

Energy balance of the core-envelope. — Since the variations 
with <7p seen for the remnant as a whole are almost entirely 
due to variations in the core and envelope rather than in the 
disk, the trends we find for the core-envelope are very similar 
to those we found above for the entire remnant. 



= 0.28^p (±0.02), 



"^deg 



:0.94-0.32^p 



(±0.02). 



1/2 

Note the dependency on qp, rather than on qp as was found 
for the entire remnant. See Fig.[9]b). 

Energy balance of the disk. — For the disk, we find very little 
dependence on qp, consistent with the idea that most of the 
changes in the partitioning of energy have to do with increased 
mixing between the donor and accretor, which affects the core 
and envelope much more than the disk. Averaged over all 
mergers, we find 



(10) 
(11) 



^disk 



=0.74 



7 disk 



^disk 



= 0.19 



17 disk 

4t = 0.07 

^disk 



(±0.03), 
(±0.02), 
(±0.02). 



(12) 
(13) 
(14) 



1.0 



0.6 



0.4 



0.2 



0.0 
1.0 




0.6 

0.4 

0.2 

0.0 
1.0 



0.6 



0.4 



0.2 



0.0 



b) 



I := s- 



c) 



• • • • • 



0.0 



0.2 



0.4 0.6 

Pc.d/Pc,a 



1.0 



Figure 9. Partition of energies in (a) the overall merger remnant, (b) the 
remnant core plus envelope, and (c) the remnant disk. In each panel, the 
fraction of total energy carried in degeneracy (triangles), thermal (squares), 
and rotational (circles) energy is shown. Colors represent different accretor 
masses, encoded as in Fig.R] Triangles represent the 1.0 - 1.0A/q merger. 

Hence, the disk is composed of non-degenerate, primarily 
rotationally-supported material. See Fig. |9]c). (Note that we 
do not try to define a ratio of internal to potential energy of 
the disk or core-envelope, since the potential energy of either 
is not straightforward to determine.) 

3.2.5. Temperature and Thermal Energy 

Since heating of the remnant is achieved through shocks 
and viscous dissipation, the most heavily mixed regions 
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should also be the hottest. We focus on equatorial thermody- 
namic values, but consider the z-axis as well for similar-mass 
mergers. 

The maximum temperature. — We find that the maximum tem- 
perature on the equatorial plane, Tmax, scales with the potential 
of the accretor (Fig.[8]e).), 



GM^mjR, 



= 0.20 (±0.03). 



(15) 



This scaling is natural in the limit of highly dissimilar-mass 
merger - for each nucleon, of order GM^m^/R.^ is liberated 
and converted into thermal energy. The temperature and ther- 
mal energy profiles in Fig.[3]show that with increasing qp, ad- 
ditional thermal energy is deposited into the remnant, but this 
energy is spread over a larger region, such that the maximum 
temperature remains roughly the same even as qp approaches 
unity. 

For a dissimilar-mass merger, the highest temperature along 
the rotational axis, T^.^^, is found at the tenuous outer edge of 
the hot envelope. It is slightly higher than the maximum tem- 
perature found in the equatorial plane. With increasing qp, 
however, the difference increases noticeably due to the two 
off-center hot spots found along the rotational axis in similar- 
mass mergers. Fitting T^^^, we find kT^^^/{GM.jnj,/R.^) = 
0.24 ±0.03, though this does not capture the upturn for simi- 
lar masses well. 

All remnants with qp > 0.8 have convectively unstable 
cores along the equatorial plane. Artificially mixing these 
cores to make them isentropic decreases their maximum tem- 
peratures by 10 - 50% (not shown in Fig.|8] but see Fig.[T7]i. 
All remnants are stable against convection along the rotational 



The density at the point of maximum temperature. — For 
dissimilar-mass mergers, the density at the hottest equatorial 
point, p{T^ax), depends mostly on the donor (see Fig.[8]f).). 
For qp < 0.5, we find 



Pc,d 



= 0.34 (±0.02). 



(16) 



This proportionality again suggests that, at least for 
dissimilar-mass mergers, the donor mixes with the accretor 
up to a fraction of the central density of the donor, as alluded 
to earlier 

At <7p > 0.6, the dependence becomes less obvious, with 
pijtnax) varying from ^25 - 90% of p^. Since for these den- 
sity ratios, the donor material starts to penetrate the central 
core of the accretor - and the accretor starts to disrupt as well 
- the simple picture of the donor mixing up to a fraction of 
its own central density may be breaking down. Furthermore, 
part of the spread in density reflects that for high qp the equa- 
torial temperature profiles become nearly flat (see Fig.[3j left 
column), thus increasing the sensitivity to noise in the deter- 
mination of the location (but not the value) of maximum tem- 
perature. This also affects our results for the enclosed mass, 

Menc(7;T,ax), bcloW. 

Since for dissimilar-mass mergers, T^^^ is located near the 
tenuous outermost regions of the hot envelope, where particle 
noise is high, the density at the point of maximum rotational 
axis temperature, p(T^.^-^) (plus symbols in Fig. [8]f).), varies 
wildly between about p(T^ix) and one order of magnitude be- 
low piT^ax)- For similar-mass mergers, p(T^^^) appears to be 
~20-50%of/o(7;n„). 



The mass enclosed within the radius of maximum temperature. — 
For dissimilar-mass mergers, the radius of maximum temper- 
ature occurs at an enclosed mass of Menc(7Inax) — A/a, while 
for mergers with ^ 0.8, maximum temperature occurs at 
the center and Menc(7max) — (see Fig. |8]g).). For qp < 0.5, 
we find 



= l-0.28^p (±0.01), 



(17) 



where the fit's y-intercept is forced to unity. Note that 
maximum temperature occurs near the bottom of the mixed 
zone, which is why Menc(7inax) is substantially smaller than 
Mi^nc{\Mi). The reasons it starts to deviate from a tight trend 
at qp ~ 0.6 are the same as those for piJ^ixY the break-down 
of the simple mixing picture and the difficulty in determining 
the location of peak temperature for a broader temperature 
plateau. 

Since Mi.^^(T^^^^) < M.^, it may be surprising that piT^^x) is 
not higher than p^.d- This is because the additional thermal 
and rotational support against gravity reduces the density gra- 
dient that would be required if degeneracy pressure were the 
only source of support. 

For dissimilar-mass mergers, Menc(7"max) is only slightly 
higher than Menc(rmax), since the remnant core is nearly spher- 
ically symmetric, apart from the fact that the hot envelope is 
slightly more extended in the vertical direction. For similar- 
mass mergers, the difference increases, reflecting the develop- 
ment of the off-center hot spots, until Menc(rjJax)/^a — 0-25 
for equal-mass mergers. 

The mass enclosing half the remnant thermal energy. — As a more 
robust measure of where thermal energy is deposited during 
the merger, we consider the mass enclosing half the remnant 
thermal energy, Menc(^£th) (see Fig. [sjh).). We find this is 
very close to the mass of the accretor. 



enc 



= 1.06 (±0.04), 



(18) 



if the 1.0- 1.0 Mq merger is neglected. One sees turnovers at 
the extremes, for qp < 0.2 and qp > 0.8. The former likely is 
because thermal energy is deposited into a narrow strip right 
on the surface of the accretor, while the latter is probably due 
to the disruption of the accretor. 

While smaller than M^^^i l^M^), Menc(^£'th) is always larger 
than Menc(7lTiax)- This reflects that high density degenerate 
material has lower specific heat, so that for the same energy 
per unit mass the temperature is higher (see Fig.|3]l. 

The width of the remnant thermal energy. — The mass enclosed 
between the 25* and 75* percentiles of thermal energy, 

AMenc(£'th) = A^enc(5£'th)-A^enc( j£'th), is a mcasurc of the ex- 
tent of the remnant that has been heated (see Fig.[8]i).). Ignor- 
ing the 1.0- 1.0 Mq merger, it can be fit by. 



AMenc(£th) 
Ma 



:0.11+0.94^p 



(±0.03). 



(19) 



Here, we did not force the y-intercept to go to zero, which is 
expected physically but gives a substantially poorer trend. 

3.2.6. Angular Velocity and Rotational Energy 

For a dissimilar-mass merger, the donor carries most of the 
angular momentum. As a result, the hot envelope contains 
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more angular momentum and features higher angular veloc- 
ities than the core, since the envelope is where most of the 
accreted donor material resides. Spin-up of the accretor is ac- 
complished through shocks and viscosity. For a similar-mass 
merger, the two stars carry similar amounts of angular mo- 
mentum and thoroughly mix. Conservation of angular mo- 
mentum then implies that the entire remnant rotates rapidly. 

The maximum angular velocity. — On the equatorial plane, the 
highest angular velocity, fimax, scales linearly with the orbital 
angular velocity of the pre-merger binary, fJoib = 2tt //"orb (see 
Fig. §]).), 

o 

= 3.8 (±0.6). (20) 

i 'orb 

The ratio of core-envelope to total angular momentum. — For 
more similar-mass mergers, more angular momentum is de- 
posited in the accretor and ends up in the core and envelope 
(see Fig. [8]k).). The ratio of core-envelope to total angular 
momentum, Lce/^^tot, is approximately. 



= 0.70^p (±0.03), 



(21) 



where we fit only for qp > 0.25 and ignore the 1.0 - 1.0 Mq 
merger For qp < 0.25, the trend becomes shallower, result- 
ing in a non-zero intercept. This suggests that even in cases 
where the donor has negligible mass some angular momen- 
tum is transferred to the accretor 

The mass enclosed inside the radius of maximum angular ve- 
locity. — For dissimilar-mass mergers, both Menc(rimax) and 
Menc(7'max) are about equal to Ma, with Menc(f^max) sUghtly 
larger than Menc(7max): Me„c(l^max)/Ma = 1.05 ±0.06. This is 
consistent with the idea that the hottest and most spun-up re- 
gions are those where the donor mixed most strongly with the 
accretor For qp ~ 0.6, the off-center angular velocity peak 
is replaced by a plateau, and Menc(f^max) becomes ill-defined; 
for even larger qp, the highest velocities occur in the center, 
and Me„c(f^max) ^ 0. See ¥'\g.M 



The mass enclosing half the remnant rotational energy. — Like for 
the thermal energy, for very dissimilar masses, the mass en- 
closing half the rotational energy, Menc(5£rot), is similar to the 
accretor mass (see Fig.jsjm).). For qp < 0.8, we find 

^enc(o^rot) 



Ma 



1.12 + 0.27^'/^ 



(±0.01) (22) 



Note that unlike Menc(i^max), A/enc( j^rot) continues to increase 
with qp (except for exactly equal-mass mergers), a conse- 
quence of particles with lower angular velocity but large lever 
arm that carry substantial rotational energy (see Fig.[3]l. Near 
qp ~ 0.8, the trend breaks as both stars are significantly dis- 
rupted. However, even exactly equal-mass mergers have more 
of their rotational energy stored in the outskirts (otherwise one 
would have M^n^{\E,odlM.^, ~ 1). 

The width of the remnant rotational energy. — We measure the 
extent to which the remnant is affected by spin-up through 
the difference between the masses enclosing 25 and 75% of 

the rotational energy, AMenc(£'rot) = Me„c( j£'rot)-Menc(3£'rot) 
(see Fig.lSln).). Ignoring the 1.0 - 1.0 Mq merger, it is well- 
describedby. 



Like for the thermal energy, one sees that for more similar- 
mass mergers, rotational energy is spread more widely 
throughout the remnant. 

3.3. A Qualitative Picture of the Merger 

From our empirical results above, a qualitative picture of 
a merger emerges. A dissimilar-mass merger has the donor 
overflowing its Roche lobe and forming an accretion stream. 
This stream mixes with the accretor up to approximately the 
central density of the donor pre-merger, pc d- Those layers of 
the accretor that are denser than p^.d are hardly affected, and 
form the cold core of the merger remnant, while the mixed 
material will form a partly-thermally supported outer enve- 
lope, which somewhat compresses the core, as well as a rota- 
tionally supported disk. 

At qp > 0.6, the above picture begins to break down, as 
portions of the donor start to penetrate to the center of the 
accretor This results in substantial heating and spin-up of 
the central core: the merger becomes a similar-mass merger. 
As the masses become more similar, the distinction between 
donor and accretor is lost and both stars disrupt and form ac- 
cretion streams. For all qp > 0.6, the remnants are similar: a 
large, ellipsoidal and partly rotationally supported warm core 
with two hotspots off the equatorial plane, surrounded by a 
small, hot disk. 

For all mergers, the maximum temperature reached by dis- 
sipation of orbital energy is proportional to the accretor' s 
gravitational potential energy. For increasing qp, the maxi- 
mum temperature remains similar, but the region over which 
the thermal energy is deposited widens. The density at max- 
imum temperature is of the same order of magnitude as the 
central density of the donor, consistent with the mixing pic- 
ture discussed above. The latter no longer holds for qp > 0.6, 
when the entire remnant is mixed and heated. 

For a dissimilar-mass merger, the angular momentum re- 
mains in the outer regions, since most of it was originally 
carried by the donor Angular momentum can be transferred 
between regions by shocks or viscosity, which becomes in- 
creasingly important as the donor penetrates deeper. As a 
result, with increasing qp, the remnant core is spun up fur- 
ther Where both WDs disrupt, leading to colliding accretion 
streams, even the densest regions of the remnant have high 
rotational velocities. 

4. VARIATION OF MERGER PARAMETERS AND ROBUSTNESS OF 

RESULTS 

In our parameter space study, we focused on the effects of 
varying the masses of the two WDs, fixing the initial separa- 
tion flo, merger completion time criterion, and WD composi- 
tion. To determine how robust our results are, we ran simula- 
tions varying these assumptions. 

4.1. Changing the Composition 

Ignoring fusion, WD mergers should be insensitive to 
changes in composition, since the dominant electron degener- 
acy pressure only depends on the mean molecular weight per 
electron, which is close to /if ~ 2 for all likely compositions. 
To confirm this, we ran simulations assuming pure helium and 
pure magnesium for an equal-mass case (0.4 - 0.4 M0) and an 
dissimilar one (0.4 - 0.8 Mq). The results are shown in Figure 



AMenc(£rot) 
M,, 



= 0.12 + 0.77^p (±0.03). 



(23) 



10 One sees that most quantities indeed have very similar 
profiles. 

The set of profiles showing most variation are those of the 
temperature. These have similar shape, but different normal- 



15 



0.4 - 0.4 M,, 



0.4 - 0.8 M,- 



10 

lo' 

10^ 

10* 

10^ 

0.8 

0.6 

0.4 

0.2 

0.0 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



10 



to 10" 

to 
u 
m 



10 
10^ 



10" 



10 







— 1 — \ — 1 — \ — - 


' \ ^ 
> \ 






~ — ^ 




— 1 — \ — 1 — \ — 
^^^^^ 


- 

' J\ - 






— 1 — \ — 1 — \ — 

Si- "•-■■.V-^- 


— 1 — 1 — 1 — 1 — 










1 \ 1 \ 


/I 


^^_====^ 








1 \ 1 \ 









0.0 0.2 



0.4 0.6 



0.8 0.0 0.2 



0.4 0.6 



1.0 



Figure 10. As Fig.|3] but for 0.4 - 0.4 M© (left) and 0.4 - 0.8 Mq (right) 
mergers with different compositions: pure '^He (red), CO (orange), and pure 
^^Mg (lime). Dash-dotted lines represent profiles along the rotational axis 
rather than the equatorial plane. 

ization. Since the thermal energy curves are very similar, it 
is clear that this reflects differences in heat capacity, which 
does depend on composition: for He composition, there are 
more (non-degenerate) ions than for our standard CO mix- 
ture, boosting the heat capacity (and thus lowering the tem- 
perature for given thermal energy), while for Mg composi- 
tion, there are fewer, lowering the heat capacity (and increas- 
ing the temperature). As a result, the maximum equatorial 
plane temperatures for the 0.4 - 0.4 simulations are 0.95, 
1.48 and 1.68 x 10*^ K for He, CO, and Mg, respectively, 
while for the 0.4 - 0.8 M© simulations, they are 2.92, 4.05, 



and 4.47 X 10^ K. 

The smaller differences seen for the other profiles reflect 
small differences in initial conditions. All WDs are con- 
structed assuming T = 5 x 10*' K throughout, which implies 
more thermal energy for higher heat capacity. As a result, 
the relaxed He and Mg WDs are slightly larger and smaller, 
respectively, than the CO WD. These slight differences in ra- 
dius translate into differences in initial separation, which in 
turn cause small differences in the angular velocity and rota- 
tional energy curves. 

4.2. Varying the Initial Binary Separation 

For our simulations, we chose an initial orbital separation 
ao for which a co-rotating donor would fill its Roche lobe. 
Since our (non-rotating) WDs are equilibrated in isolation, 
once the simulation starts they immediately begin to adjust to 
the tides and hence disrupt quickly. Ideally, one would allow 
them to adjust to the binary potential and start mass transfer 
properly. For non-synchronous rotation, however, this is not 
straightforward (see Sec. \5.2\ . Nevertheless, we try to get a 
sense of the influence of this by running simulations for two 
cases - 0.6 - 0.6 Mq and 0.6 - 0.8 Mq - with aq increased and 



decreased by 10% (see Fig. 11 



Our default simulations were considered complete at 6 or- 
bits of the initial binary. For runs where aq was changed, we 
used 2.5% non-axisymmetry, and a requirement for the den- 
sity to be highest at the remnant's center, as the completion 
critericj^ Not surprisingly, runs with larger aq needed longer 
to achieve these criteria: with a 10% increase, the 0.6 - 0.8 
Mq merger required 861 seconds, or ^ 15 orbits, to complete, 
while the one with a 10% decrease required 230 s, or ^5.5 or- 
bits. Similarly, the 0.6 - 0.6 Mq merger with a 10% increase 
in Ao required ^+13 orbits (725 s) to complete, compared to 
the 7 orbits of a 10% decrease (283 s). In both cases, the in- 
crease simply reflects that it takes longer for the donor to be 
disrupted fully if aq is increasecQ For instance, for the 0.6 
- 0.8 Mq binary with increased separation, it took almost a 
dozen orbits before full disruption, while at the standard sep- 
aration disruption occurred after just 1.5 orbits. This feature 
is of particul ar interest because for m erg ers of synchronously 
rotating WDs |Dan et"aL] ( |20TT][20T2l ) and |Raskin etaL|pOT2) 
all note almost immediate disruption of the donor when ap- 
proximate initial conditions are used, and much delayed dis- 
ruption for m ore accurate initial conditions (up to ^30 orbits; 
see Sec. |5.2| i. 

We find that the density profiles of the merger remnants are 
remarkably insensitive to varying aq (pc changing by < 2% 
for the 0.6 - 0.8 Mq merger, and - 20% for the 0.6 - 0.6 M© 
merger), and show substantial systematic changes only in the 
outer regions. The latter can be understood from the mixing 
and rotational profiles, where one sees that with increasing aq, 
donor material is mixed less deeply into the accretor, and ro- 
tational energy is shifted outward, causing the rotational fre- 
quency to peak at lower values and larger radii. This reflects 
the increase in angular momentum with increasing aq, which 
creates a more rotationally supported remnant. In the 0.6 - 
0.8 M0 merger, the decreased mixing causes the accretor to 
be spun up less, thus lowering the rotational energy of the 

^ The 2.5% non-axisymmetry convergence time is 312 s (6.6 orbits) for 
our default 0.6 - 0.6 Mq , and 250 s (5.2 orbits) for our default 0.6 - 0.8 Mq . 

' Of course, if placed far enough, the binary does not merge. For a 0.6 - 
0.8 Mq binary, no mass transfer occurred within 500 s if ao was increased by 
20%. 
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Figure 11. As Fig.[T()] but for 0.6 - 0.6 M© (left) and 0.6 - 0.8 Mq (right) 
mergers with varying initial orbital separation: 0.9 (red), 1.0 (green), and 1.1 
(blue) times the value used for the parameter space study. 



Figure 12. As Fig.[To] but for 0.6 - 0.6 M© (left) and 0.6 - 0.8 Mq (right) 
mergers, comparing our default, irrotational case (blue) with that assuming 
synchronous rotation (red). 



core, and narrowing the thermal energy plateau. These effects 
are also seen in the similar-mass case, where the center of the 
remnant receives less rotational support and becomes denser 
with increasing separation, and the mixing becomes less uni- 
form. 

Qualitatively, with increasing aq, the properties change in a 
way that is similar to the changes seen with decreasing qp, i.e., 
similar to mergers with more dissimilar mass: reduced mix- 
ing, larger disks and less core rotational support, and shifts in 
the thermal and rotational energies toward larger radii. The 
converse is also true, decreasing aq has similar effects as in- 
creasing qp, i.e., the mergers become similar to those with 



more equal masses. The changes are substantial at times: e.g. 
with a 10% increase in aq for the 0.6 - 0.6 Mq remnant, the 
maximum equatorial temperature is reduced by 40%, while 
the corresponding density increases by 25% (for the rotational 
axis hotspots, the values are a 13% and 30% reduction, re- 
spectively), and the mass of the disk increases by 65%. Sim- 
ilar, though less extreme, changes are seen for the properties 
of the 0.6 - 0.8 Mq remnant. All this makes ao one of the 
parameters our mergers are most sensitive to. 

4.3. Synchronization 
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In our simulations, the WDs have zero spin, i.e., we assume 
that tidal dissipation is too weak to synchronize their rotation. 
Whether or not this is correct is currently unknown, but to see 
what the effect could be, we ran simulations assuming syn- 
chronized rotation for 0.6 - 0.6 Mq and 0.6 - 0.8 Mq binaries 
(see Fig. 



0.6 - 0.6 M 



0.6 - 0.8 M,- 



12i 



As in Sec. 4.2 our unsynchronized runs are from our pa- 



rameter space study, and use 6 orbits as their completion crite- 
rion, while our synchronized runs use 2.5% non-axisymmetry, 
and a requirement for the density to be highest at the rem- 
nant's center Completion occurred at 407 s (^8.5 orbits) for 
the synchronized 0.6 - 0.6 Mq simulation, and 314 s (^6.5) 
orbits for the synchronized 0.6 - 0.8 Mq. 

The asynchronous and synchronous mergers differ mostly 
in the amount of heating and spin-up. This has two causes. 
First, for the synchronized binary, the total amount of angular 
momentum is about 10% larger, and high angular momentum 
material has more difficulty penetrating the accretor, as is ev- 
ident in the 0.6 - 0.8 Mq mixing profile (because of the larger 
amount of angular momentum, the mergers also take about 
1.5 orbits longer to achieve 2.5% non-axisymmetry). Second, 
in a synchronized binary, the donor and accretor have much 
less differential rotation with respect to each other, leading to 
much less spin-up and heating due to shocks and shearing. 

Both effects are largest for the equal-mass case. In partic- 
ular, in a synchronized, equal-mass binary contact can occur 
without any friction, while in an unsynchronized one it in- 
volves shocks at the full orbital velocity. In consequence, for 
the synchronized case, rotational support is weaker in the cen- 
ter and stronger in the outskirts, causing the central density 
of our 0.6 - 0.6 Mq remnant to increase by ^ 70% and the 
disk mass to increase by a factor of 2. Furthermore, while for 
the non-synchronized case, the maximum temperature along 
the equatorial plane was found in the center, for the synchro- 
nized case it is found in the outskirts, and is more than a factor 
two lower (1.3 x 10** K instead of 2.9 x 10** K). The hotspots 
on the rotational axis also have much reduced temperature, 
2.3 X lO** K instead of 3.6 x 10^ K. 

For the dissimilar-mass merger, the effects of synchroniza- 
tion are less dramatic: the accretor still spins up substantially, 
and rotational and thermal energy are deposited in roughly the 
same way. The main difference is that the synchronized case 
has slightly less mixing, causing a drop in total thermal energy 
and maximum temperature (from 4.1 x 10** K to 3.5 x 10** K 
on the equatorial plane, and from 4.6 x 10'^ K to 4.3 x 10^ K 
on the rotational axis). 

4.4. Running the Simulation Longer 

We considered our mergers completed after 6 orbits, since 
in that time they on average had reached our convergence cri- 
terion of 2.5% non-axisymmetry (see Sec. |2.4[ ). To test the 
robustness of our results, we also determined properties at- 
tained after 8 orbits. In Fig. [13] we compare the 6 and 8 orbit 
results for 0.6 - 0.6 Mq and (l6 - 0.8 M© binaries. 

We find our mergers show little evolution between 6 and 
8 orbits, with the largest changes seen for the angular veloc- 
ity profiles. For the 0.6 - 0.8 Mq merger, the rigidly rotat- 
ing core sped up and the off-center peak decreased in height 
and moved out, while for the 0.6 - 0.6 Mq merger, the cen- 
ter spun down and the rotational profile became flatter. In the 
dissimilar-mass merger, Tmax and /o(7'niax) changed by < 5% 
and the density and temperature structures nearly overlap, 
while in the equal-mass merger Tmax and p(T^^) changed by 




0.4 0.6 



0.4 0.6 



Figure 13. As Fig.[To] but for 0.6 - 0.6 Mq (left) and 0.6 - 0.8 Mq (right) 
mergers, comparing properties for our default simulation time of 6 initial 
orbital periods (blue) with those obtained after 8 orbital periods (red). 

-20% (2.9 to 2.3 X lO'^K and 1.7 to 2.0 x lO^gcm'^), re- 
flecting an increase in central density and a shifting of the 
temperature profile, with temperature decreasing in the cen- 
ter but increasing elsewhere. The evolution of all properties 
is consistent with viscous evolution - expected to follow the 
merger proper - with the core driven into rigid rotation, and 
angular momentum transferred outward to the disk. In the 
dissimilar-mass merger, the net effect is spin-up of the core 
and spin-down of the envelope, while in the equal-mass case 
it is the reverse. Of course, in the process, rotational energy is 
turned into thermal energy, heating the remnants. 
One curious aspect for equal-mass mergers is the evolution 
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of the off-center hot spots (Fig.|2|i. Over time, these broaden 
parallel to the equatorial plane, yet become narrower along 
the rotational axis. As a result, the hourglass shape is lost, 
and the center of the remnant stops being one of the hottest 
point in the system. After 8 orbits, the system resembles more 
closely what we find for typical similar-mass mergers, which 
have more pancake-shaped hot spots flanking a colder, denser 
region on the equatorial plane. 

Comparing more broadly the 6 and 8-orbit results, we find 
changes at the ^5% level. The trends presented in Sec. |3.2| 
continue to hold to within ^ 10% except for (which be- 
comes ^ = 0.98 + 0.74^2 (±o.07)), the fraction of disk en- 
ergy in degeneracy energy (£'deg,disk/£'i,di,sk = 0.07 (±0.02)), 
mass enclosed within the radius of maximum temperature 
(Menc(7lTiax)/-^a = l-0.21^p (±0.01)), maximum rotational 
frequency (fimax/f^orb = 3.4 (±0.5)) and the widths of the re- 
gions in which thermal and rotational energy are deposited 
(AMe„c(£th)/Ma = 0.87^^ + 0.13 (±0.03); ^M, I M, = 
0.70^p+0.15 (±0.02)). Changes to these trends are consistent 
with the viscous evolution described above: the remnant is be- 
ginning to spin down, lose its rotational support, and energy 
is being redistributed. 



4.5. Viscosity Prescription 

The addition of artificial viscosity is required in SPH to ac- 
curately capture shocks, but no consensus exists on how best 
to implement it. We ran two additional simulations of a 0.6 
- 0.8 M0 merger to check the robustness of our results with 
respect to changes in the viscosity, one with small and one 
with large artificial viscosity (fixed (a,/?) = (0.05,0.1) and 
(1,2), respectively). Here, we expect that low values of a 
will lead to large particle noise and inaccurate shock captur- 
ing, while high values result in large viscous heating and rapid 
loss of differential rotation. Our results confirm this (Fig. 14 1: 
the simulation with low artificial viscosity leads to a remnant 
with stronger differential rotation, with the disk carrying 34% 
more rotational energy (and the remnant 33% less) than in the 
standard variable a simulation. Lower viscosity also leads to 
greater mixing of donor and accretor material, reflecting the 
stronger diffusion associated with the larger particle noise in- 
herent to low viscosity. 

Aside from the mixing, the results for the three different 
viscosity prescriptions do not differ greatly. While one might 
have expected greater dissipation of rotational into thermal 
energy for higher viscosity, the maximum temperatures and 
rotation rates vary by < 10%, and the thermal and rotational 
profiles are quite similar. The density profiles are virtually 
identical except near the outer parts, where the low viscosity 
simulation leaves matter with greater rotational support. 



Viscosity Test Numerical Resolution Test 



4.6. Spurious Heating 

As discussed in Sec. 12.1 1 and throughout Sec. |3] SPH 
particle-particle interactions lead to small increases m inter- 
nal energy, which, in the most degenerate regions, correspond 
to significant temperature increases. Here, we discuss the ex- 
tent to which this "spurious heating" affects our results. 

As a comparison for the spurious heating seen in some of 
the simulations, we relaxed a 0.8 isolated white dwarf for 
489 s longer than the standard 8 1 s we used for relaxing single 
stars. While the total energy of the WD (potential, degeneracy 
and thermal energy combined) increased by '-^ 1% of the orig- 
inal total energy over the extended period of time, the change 
in thermal energy was enough to raise the central WD temper- 
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Figure 14. As Fig. |10| but comparing simulations of a 0.6 - 0.8 Mq merger 
with different viscosity prescriptions (left) and simulations of a 0.625 - 
0.65 Mq merger with different numerical resolutions (right). For the vis- 
cosity, we compare fixed low viscosity (a = 0.05, /3 = 0.1; blue), standard 
variable viscosity (green), and fixed high viscosity (a = 1.0, /3 = 2.0; red). 
For the numerical resolution test, we show simulations at one quarter (red), 
half (orange), and double (blue) the default resolution, as well as the default 
simulation (lime) and a rerun of the default simulation (green) to determine 
the effect of order-of-execution differences, round-off errors and other such 
numerical effects. 



ature from 5 x 10'' K to 1.2 x 10** K (even within the standai'd 
81s, the temperature rises to 1.2 x 10^ 
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K). In Fig. 

compare the thermal energy profile of this isolated Mq 
WD with those of a 0.4 - 0.8 Mq and a 0.7 - 0.8 M© merger, 
with the times for the isolated WD taken at 489 and 224 s (the 
mergers' respective completion times) longer than the stan- 
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Figure 15. Specific thermal energy as a function of enclosed mass for a 0.4 - 
0.8 Mq (blue) and a 0.7 - 0.8 Mq merger (red), shown both along the equa- 
torial plane (solid curves) and along the rotational axis (dot-dashed). Also 
shown are the (spherical) profiles found for an isolated 0.8 Mq white dwarf 
(dashed) simulated using the same parameters, and for the same completion 
times (six initial orbital periods, equivalent to 489 s and 224 s). Spurious 
heating is responsible for nearly all the thermal energy in the core of the 0.4 
- 0.8A/q remnant, and for about one third in the core of the 0.7 - O.8M0 
remnant. It is not important in regions heated by interaction. (The "hook" in 
the outer layers of the white dwarf profile reflects the high initial temperature 
chosen; in a merger, this is erased by the interaction.) 

dard 81s. The total thermal energy generated in the WD at an 
additional 489 s is ^10% of the thermal energy generated in 
a 0.4 - 0.8 M0 merger]^ Indeed, given how well the specific 
thermal energy profiles match in the interior, it is clear that 
spurious heating dominates there. 

Spurious heating is much less important for the 0.7 - 0.8 
merger, since its core has mixed to much greater extent and 
less time was needed for the merger to complete. The total 
thermal energy generated in the isolated WD at 224 s is only 
^3% of the thermal energy generated during merger, and even 
at the very center spurious heating contributes only ^ 35% 
(rather than nearly all) of the thermal energy. As a result, the 
central temperature of the merger, 1 .4 x 10** K, is substantially 
higher than that of the isolated WD, 7.8 x 10^ K. 

Overall, we conclude that spurious heating is present, but 
is recognized fairly easily and does not influence our con- 
clusions. In particular, its effects should be small in both 
high-density regions with T > 3 x 10^ K and in low-density, 
< 10^ gem"-' material where degeneracy has been lifted. 

The spurious heating may be minimized by an appropri- 
ate choice of the viscosity prescription. For instance, a larger 
viscosity would reduce particle noise, which in turn would re- 
duce spurious heating. However, a large viscosity would also 
increase dissipation of differential rotation. Other simulations 
may suffer from spurious heating as well. In this respect, it 
is intriguing that our equatorial temperatu re curve s for a 0.6 - 
0.8 Mq merger are a good match those of LIG09| even in the 
center (Fig. 16 Sec. 5.1 1. 



A ~10% increase in thermal energy corresponds to ~ 1% increase in the 
overall energy of the 0.4 - 0.8 Mq remnant, somewhat larger than the typical 
~ 0.3% level at which Gasoline conserves total energy in our simulations. 
We find that siinilar-mass simulations tend to lose total energy at the 0.05% 
level, while some of the low qp mergers gain more than 1 % in total energy 
due to spurious heating. 



4.7. Resolution 

To determine whether or not particle resolution matters for 
our results, we ran three additional simulations of a 0.625 - 
0.65 Mq merger, with roughly a quarter, half, and double 
the resolution (63,736, 127,525 and 510,047 particles, respec- 
tively), as well as a second simulation with the same resolu- 
tion (255,035 particles). Here, we chose 0.625 - 0.65 M© to 
see if numerical resolution has any effect on whether a merger 
is "similar-mass". 
From Fig. 14 one sees that the two runs using the same 
and identical initial conditions and the 
still give slightly different results. 



number of particles 
same version of Gasoline - 
This is due to the inherent non-linear nature of fluid dynam- 
ics, coupled with small, random perturbations, e.g., from dif- 
ferences in the order of force addition in parallel processing, 
round-off errors and slight inconsistencies in converting ther- 
mal energy to temperature. Overall, merger remnant prop- 
erties change by ~ 3.5% between the two runs. The most 
prominent differences are seen in properties determined from 
low numbers of particles, such as % (varies by ^ 10%), and 
those involving finding maxima of temperature plateaus, such 
as p{T^.j,y) (40%) and piT^.^,,) (a factor of 3 - in one case con- 
vection shifts the temperature maximum off-center). 

The differences for different resolutions are larger While to 
first order, the equatorial density and mixing profiles are very 
similar, there is a systematic ^ 20% drop in the equatorial 
density - ^-^25% in the rotational axis density - near the cen- 
ter of the remnant with increasing particle resolution (top right 
panel of Fig. 14 1. The angular velocity and rotational energy 
profiles are agam very similar, except in the central regions, 
where there is a ^20% increase in ilmax- For the tempera- 
ture the effects are larger: with higher numerical resolutions, 
most of the equatorial plane is colder, with a ^ 20% drop in 
the value of the temperature plateau nearM/Mtot = 0.5. The 
temperatures along the rotational axis, however, increase with 
increasing numerical resolution, by ~ 10% across the range of 
resolutions, as does the upturn in equatorial temperature near 
the center of the remnant, by ^ 50% (^ 30% if we do not in- 
clude the lowest resolution run). The latter effects are due to 
increasing prominence of the off-center hotspots at higher nu- 
merical resolutions, which also tend to look more hourglass- 
shaped. Indeed, for our lowest resolution, the densest mate- 
rial in the two stars remains relatively cold throughout the en- 
tire merger, resembling the synchronized systems described 
in Sec. |4.3| Finally, we find that if we do not consider the 
lowest resolution run, the disk half-mass radius varies by 4%, 
angular velocity at the half-mass radius varies by 4%, and the 
core-envelope mass changes b y 3%. This is simil ar to the re- 
sults of the resolution tests of Raski n et al.| ( |2012] ). 

The 2.5% non-axisymmetry convergence times for the half 
and double-resolution runs are within 14 s of the 275 s non- 
axisymmetry convergence time of the default run, a small dif- 
ference that implies a negligible amount of post-merger evo- 
lution. Only the quarter-resolution run deviated substantially, 
converging 57 s earlier. This may simply reflect the smaller 
number of particles in the disk, where the system is most 
asymmetric. 

We stress that even though the order-of-magnitude change 
in resolution generates 10 - 30% variations in some remnant 
properties, the overall shapes of the profiles in Fig. 



14 



are 



very similar. In particular, the merger remnant does noTTook 
more or less "similar-mass" (except, arguably, the tempera- 
ture curve at the lowest resolution). We thus conclude that 



20 



our trends should be robust. 



5. COMPARISON WITH OTHERS 

5.1. Comparison With \Loren-Aguilar et a/.| ( [2009p 



LIG09| simulated a number of WD mergers, and gave de- 
tailed temperature, surface density, and rotational frequency 



curves for three. In Fig. 16 we compare their results (from 
their Figs. 3 and 4) with ours for two of these, 0.6 - 0.6 Mq 
and 0.6 - 0.8 Mq (the third was a 0.4 - 0.8 Mq He - CO WD 
merger, whose temperature profile cannot be compared di- 
rectly). We note that they used different initial conditions, 
starting their systems with an orbital separation too large for 
mass transfer to begin, and then slowly reducing the separa- 
tion until it does. This point defines their t = and the start of 
the merger simulation proper. Given this different setup, their 
merger completion times cannot be compared directly to ours. 
In their simulations, however, coalescence (the final consoli- 
dation of the two WDs into one) also occurs after just about 
one orbit, so the differences should not be too large. To give 
a sense of the effect of different completion criteria, we com- 
pare their results both with our standard results, taken after 6 
orbits, and our results taken at their merger completion times. 

For both mergers, the surface density curves are similar, 
although in their 0.6 - 0.6 M© merger, the central peak is ^ 
30% higher (^ 10% if we use their completion time of 514 s). 
For the 0.6 - 0.8 M© merger, the temperature profiles are also 
very similar, with maxim £0differing by only -10 - 15%, and 
having nearly identical shapes. Larger differences are seen 
for the 0.6 - 0.8 M© rotational frequency profile, where the 
angular velocity peaks further out and at lower value (—0.3 s"' 
compared to our 0.45s"' - or 0.50s"' using their completion 
time of 164 s). Indeed, our entire remnant is more spun-up 
than theirs. 

For the 0.6 - 0.6 M© merger, LIG09 have a plateau in their 
angular frequency profile, with Jlmax — 0.25 s"', while our 
profile is much more peaked and reaches a much higher fre- 
quency, of 0.60s"'. By their completion time, our rotation 
curve is not as sharply peaked, but still reaches 0.44s"'. The 
temperature profiles are also much less similar: at our comple- 
tion time, our maximum temperature in the equatorial plane 
is a factor of 2 lower than theirs (factor of 3.3 at their com- 
pletion time), and even our maximum temperature along the 
rotational axis is a factor 1 .6 lower (factor 1 .9 at their comple- 
tion time). 

Finally, we can compare h ow mas s is distributed. In both 
our simulations and those of 'LIG09', negligible mass is lost, 
so only the distribution between disk and core-envelope mat- 
ters. For our 0.4 - 0.8 Mq, 0.6 - 0.6 Mq, and 0.6 - 0.8 Mq 
simulations, we infer disk masses of 0. 3 1 , 0. 1 0, and 0.40 Mq , 
respectively, which are reasonab ly close to the 0.28, 0. 10, and 
0.30Mq, respectively, Usted by |LIG09| (their Table 1), espe- 
cially considering that we hkely use a different definition of 
what is "disk". 

Overall, the primary differences between our simulations 
appear to be the amount of spin-up and heating of the equal- 
mass merger. We believe it is unlikely that this reflects differ- 
ences in initial conditions: we found much smaller changes in 
the angular velocity profile with increasing ao (see Fig. [TT]l, 



and in the simulation of ILIG09I the stars still seem to be 
quite close to spherically symmetric at the start and disrupt 
quickly (their Fig. 1), even though they were more properly 
relaxed. Instead, we believe the more likely explanation is 
that the viscosity prescription of LIG09 based in Reimann 
solvers, yields larger effective viscosity. This would explain 
both the reduction in angular velocity and increase in tem- 
perature (since viscous evolution converts rotational into ther- 
mal energy), as well as the fact that s imilar-m ass mergers are 
affected more (they mix more, and LIG09| ran their equal- 
mass merger for a very long time). If we ran our simulations 
longer and thus included further viscous evolution, the simi- 
larity with their simulations would likely be closer. 

5.2. Comparison with Others 

Simulations of WD m ergers hav e also been presente d by 
Yoon et al. (2007), Pakm or et al.ir 2010 2011, 20121, [Dan] 
,et al.. ( ,20 11,, 20 12) , and Raskin et al., ( 2012) . Unfort unately 
comparison with those results is difficult, since, unhke LIG09 
all these authors are sparse with quantitative details about 

their results. 

An exception is the 0.81 - 0.9 Mq merger simulated by Dan 
et al. (2012 1, shown in thek Fig. 1. While that simulation is 
for synchronized WDs, it is still particularly useful to com- 



pare with, since Dan et al.l show results for both approxi- 
mate and accurate initial conditions. We find that their spher- 
ically enclosed mass profile is very similar to ours, with, e.g., 
M = 0.9Mq at 4.5 x 10^ cm in both (though since spherically 
enclosed mass is a cumulative quantity, significant structural 
differences can remain hidden). Our spherically averaged 
density profile looks most similar to the profile they found 
using approximate initial conditions. Our central density, 
1.9 X lO^gcm"-*, is within — 10% of theirs, and the density 
profile remains similar up to r ~ 5 x 10*^ cm (p ~ lO^gcm"-'). 
Beyond, their profile becomes shallower while ours continues 
to decline; at r = 10^ cm, they find p ~ 3 x 10^ gem"-', while 
we find p ~ 10^ gcm"^. This may be a consequence of the ad- 
ditional angular momentum associated with synchronized ro- 
tation. With accurate initial conditions, a difference with our 
results is that the density profile becomes flat beyond 10^ cm. 

Comparing temperature profiles, we roughly reproduce 
their spherically averaged one for approximate initial con- 
ditions, including the off-center peak - their r^ax is — 20% 
low er (to be expected since their binary is synchronized; see 
Sec. lijl, but is also located at 4 x 10** cm (or Mr ~ 0.9Mq). 



However, our central temperature (2.2 x 10 K) is an order of 
magnitude higher than theirs (2 x 10^ K) - almost certainly a 
result of the much larger dissipation expected for non-rotating 
WDs - and at large distances, r>lQ^ cm our temperatures are 
systematically hotter. With accurate initial conditions they 
found an even narrower temperature peak than the one with 
approximate conditions, which thus deviates even more from 
our curve. This trend is similar to what we see when increas- 
ing fl() (Sec. |4.2| i, so it seems likely we would reproduce their 
simulations more closely if we used the same initial condi- 
tions. 

5.3. The Importance of Accurate Initial Conditions 
Many of the rec ent simulations (Dan et al. 201 1[ |2012[ 



' Maximum temperatures given in Table 1 of LIG09 refer to hot spots in 
their simulations, and are about a factor of 2 higher than the hottest points on 
their temperature curves. As we have not done hot-spot finding, we cannot 
compare with those values. 



Raskin et al. 2012 1 assume co-rotating WDs. ITiis assump 



tion IS numerically convenient, in that it is relatively straight- 
forward to start the simulation in the physically correct state: 
since in the co-rotating frame there are no flow velocities, one 



21 




Figure 16. Comparison of our results with those of'LIG09' for a 0.6 - 0.6 Mq (top) and a 0.6 - 0.8 Mq (bottom) merger. Shown are surface density, remnant 
(solid) and Keplerian (dashed) angular frequency, and temperature (left to right), with profiles from LIG09 in blue, and our equivalent ones in red and orange. 
Here, the former are for our default completion time of 6 initial orbital periods and the latter for their completion time (514 s or 10.9 orbital periods for the 0.6 - 
0.6 A/0 merger, and 164 s or 3.4 orbits for the 0.6 - 0.8 merger). 



can easily relax a simulated binary within an appropriate po- 
tential in the co-rotating frame, damping out any velocities 
resulting from an initial mismatch. 

As a result, it has been possible to st udy the onset of mass 
transfer in detail. As first pointed out by D'Souza et al.|p006] l 
from simulations using a grid code, the disruption of the donor 
is preceded by a rather long - dozens of orbits - phase of 



mass transfer Further simulations by Dan et al. (201 1 2012 1 
showed that in this initial phase a significant fraction, ~ 10% 
of the donor mass, is transferred. As a result, e.g., the disk 
is substantially colder and more extended. The remnant core 
seems more subtly affected, in that its appearance becomes 
"more dissimilar", reflecting that coalescence is between two 
WDs whose masses have become more disparate than they 
were initially. As a consequence, e.g., even for similar-mass 
binaries, the hottest point of the merger is fou nd to be well 
outside the center. Indeed, |Raskin e t al.| (j20T2| find that even 
for equal-mass binaries, the final outcome for more massive 
mergers is one where the core of one of the WDs is virtually 
undisturbed. 

At present, it is not clear how important accurate initial con- 
ditions would be for asynchronous mergers. Qualitatively, we 
expect the effects to be smaller than for synchronous merg- 
ers, for thre e reasons. First, from the analytic study of Lai 
|et al.| ( |1994| l, in which tidal and rotational distortion are ap- 
proximated by ellipsoids, co-rotating binaries always reach 
contact or Roche lobe overflow before becoming dynamically 
unstable, while irrotational binaries become dynamically un- 
stable first. While an exact treatment of the irrotational case 
found that, in fact, Roche con tact preceded dynamical insta- 
bility (|Uryu & Eriguchi|1998[ ), it suggests that WDs in irrota- 
tional binaries will disrupt much sooner Second, the simula- 



tions of LIG09 use initial conditions that should b e qui te close 
to correct, yet their WDs disrupt quickly (see Sec. |5.1|l. Third, 



the two components are counterrotating in the rotating frame. 
Hence, any mass transferred will hit the accretor with a larger 
relative velocity than would be the case for co-rotating WDs. 



Indeed, in the limit of equal-mass WDs, very little would hap- 
pen for co-rotating WDs when one reaches contact, while a 
strong shock would be expected for the irrotational case. In 
general, one expects part of the shocked material to enter a 
high- entropy halo around the accretor. For co-rotating WDs, 
|Dan et al.| \20l 1 ) found that this halo helps remove angular 
momentum from the orbit, leading to a shorter start-up phase. 
For the irrotational case, given the stronger expected shocks, 
the start-up phase would likely be reduced even further 

Ideally, one would still simulate the initial mass transfer 
phase accurately. Unfortunately , even though the equilibrium 
solution is known (|Uryu & Eriguchi 1998 ), it is not straight- 
forward to set up the initial binary properly, since it is diffi- 
cult to relax to a state that includes substantial fluid motion, 
and to slowly evolve such a state to contact, while ensuring 
viscosity remains low enough that there is no artificial tidal 
dissipation. Such dissipation is seen in our tests with varying 
ini tial d istance aq in Sec. |4.2| (and may affect the simulations 
of ILIG09 as well). Prior to coalescence, strong dissipation 
of tidal bulges heats the outer envelope of the donor (both 
stars for similar-mass mergers), and spin-orbit coupling due 
to both tides and the direct-impact accretion stream result in 
both donor and accretor becoming 25 - 50% synchronized by 
coalescence. 

Since it significantly affects the merger and merger out- 
come, whether or not tidal dissipation causes real CO WD 
binaries to synchronize before the merger remains a major 
source of uncertainty. For the radiative stellar envelopes ap- 
propriate for WDs, tidal dissipation is expected to be inef- 
ficient, with a timescale 10^^ to 10'^ yrs, suggesting that 
WDs do not synchronize ( [Marsh et al.||2004| and references 
therein). However, coupling of the tides to pulsa tions may 
dramatically increase dissipation ( Fuller & Lai|2012| ). Fortu- 
nately, it may be possible to determin e the rate o f synchro- 
nization observationally. For instance, Piro (201 1) suggested 
that tidal dissipation is responsible for the relatively high tem- 
perature of the primary WD in the 13 -minute eclipsing bi- 
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Figure 17. Merger remnant maximum temperature r,„ax and con'esponding 
density p(rmax) for all merger remnants. Values along the equatorial plane 
are marked with circles, with lines connecting points with the same accretor 
inass, while values along the rotational axis (only plotted for similar-mass 
mergers) are marked with squares (for all. colors indicate accretor mass, en- 
coded as in Fig. [5}. For similar-mass mergers, equatorial temperatures have 
been adjusted to account for mixing in coiivectively unstable cores. Also 
shown are contours of constant neutrino cooling timescale Ty = CpT jsv and 
carbon fusion heating timescale Tcc = CpT /ecc, both in years, as well as en- 
tropy S in ergK"' . (Here, Cp is the heat capacity at constant pressure and e 
the specific energy loss/gain rate.) The lines labelled Tcc = and r^c = rjy,, 
denote where the carbon fusion heating timescale balances the neutrino cool- 
ing and dynamical timescales, respectively. Finally, the P = 2P(T=Q) line is 
shown as an approximate upper bound of the region where degeneracy pres- 
sure dominates. All quantities were calculated using MESA iPaxtonetakj 
|20TT) . 

nai-y SDSS J065133.33H-284423.3, predicting that it would be 
about halfway to being synchronized. This could be tested 
by either measuring the rotational broadening of the narrow 
cores of the hydrogen lines, or looking for velocity deviations 
through the transit of the more massive secondary. 

6. POST-MERGER EVOLUTION 

We now turn to the question of how our merger remnants 
will evolve. To set the stage, we show in Fig. [iTjfor all rem- 
nants the maximum temperature T^ax found along the equato- 
rial plana^as a function of the corresponding density piT^^^^). 
Here, forme similar-mass mer gers fo r which we found con- 
vectively unstable cores (Sec. 3.2.5 i, we show the (lower) 
temperatures reached after artificially mixing them. For those 
mergers, the much higher temperatures reached along the ro- 
tational axis are shown as well (squares). One sees the trends 
identified earlier: TJ^ax is mostly set by the accretor, while 
p(7max) depends more strongly on the donor. As a result, max- 
imum temperature occurs in less degenerate conditions for 
dissimilar-mass mergers, crossing the degeneracy line for our 
most disparate cases. One also sees that for all but the most 
massive accretors, carbon fusion will not start: the neutrino 
cooling time is shorter than the fusion heating time. This is 
consistent with what was found in previous work (see Sec.[T]i. 

The central temperatures for the 0.625 - 0.65 Mq and 1.0 - 1.0 Mq 
mergers are ~4% and 10% lower than their respective maximum tempera- 
tures. In both cases, however, the center is much denser than the off-center 
hotspot, and since our estiinated post-merger evolution more greatly affects 
central material, we show the central equatorial density and temperature for 
these two systems in Fig. |17| rather than the maximum. 



Figure 18. As Fig. |17| except for maximum temperatures and corresponding 
densities following estimated post-merger evolution. The estimate assumes 
that the remnant spins down completely, that all rotational energy is used to 
drive matter to large distances, and that the remainder adjusts adiabatically 
(see text). 



6.1. Viscous Evolution and Possible Spin Down 

Following the merger, processes that happen on timescales 
slower than the dynamical time can become important. These 
include viscous evolution, neutrino emission, radiative or 
convective thermal adjustment, and magnetic dipole radia- 
tion spin-down. Out of these, convection acts on the fastest 
timescale, and we already included its effect on the core in 
Fig. [17] Next fastest would almost certainly be viscous evolu- 
tion. T^he merger re mnant is unstable to both the magneto- 
rotational instability (Balbus & Hawley"1991) and Tayler- 
Spruit dynamo (Spruit 2002). Radiative adjustment is ex- 
pected to be much slower, except at the surface, where ra- 
diative losses may also lead to convection in some system s 
dShen et al.||20T2l ISchwab et all|20T2l [Raskin et al.||20T2l ). 



Using the standard [Shakura & Sunyaevf ( |1973| l g-prescription 
for the viscosity v = aCsH, where Cj is the local soundspeed 
and H is the scale height of the system, the viscous evolution 
timescale for the remnant disk is 
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(24) 



implying a timescale fvisc ^ lO"' - 10^ s for a^lO "'-lO ' and 
fdyn ^ 10 s. This is orders of magnitude smaller than both the 
neutrino loss timescale (> lO"* y rs; see Fig. 17 1 and thermal 



Shen et al.|2()12| . 

differential rotation dur- 



adjustment timescale (> 10 yrs; 

It is possible that the strong 
ing a merger results in substantial amplification of mag- 
netic fields. The one known probable WD merger remnant, 
RE J03 17-853 has a surface magnetic field of 3.4 x 10*^0 
( [Barstow et al.|1 995; Kiilebi et al. 2010). If mergers lead to 
strongly magnetized WDs, and these WDs additionally drive 
an ionized outflow, the magnetic coupling between the out- 
flow and the WD could serve to transport angular momentum 
out of the system, spinning down the WD. The timescale for 
such a spin-down is roughly given by. 



^msd ' 



(M17)3/5(B2;;6)2/5 



(25) 
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Figure 19. Estimate of post-merger viscous evolution for 0.4 - 0.8 Mq (top) 
and 0.6 - 0.6 Mq (bottom) mergers. In blue are shown the temperature- 
density structure of the merger remnant, on the equatorial plane before (dot- 
ted) and after (solid) correction for convection, as well as along the rota- 
tional axis (dot-dashed), with points marking the hottest locations (circles and 
squares) and steps of 0.2 Mq in spherical enclosed mass (triangles pointing 
up and down). In red, estimates of the structure following viscous evolution 
are shown, where it is assumed that the remnant spins down completely, that 
all rotational energy is used to drive matter to large distances, and that the re- 
mainder adjusts adiabatically (see text). For reference , als o shown as dotted 
curves are the contours of constant entropy from Fig. [TT] (green), as well as 
the lines where Tcc = (magenta), Tcc = T^y„ (red), and y = 2P{T=0) (black). 

where L is the angular momentum of the remnant, M the 
mass loss rate, the Alfven radius, and O the angular 
spin frequency. For the second approximation, we used 
that Ra ^ (B^R'^/Mfiy^^, with B the surface magnetic field 
and R the remnant radius. Scaling to B = BglO^G and 
M = M_7lO~^MQyr~' (similar to what is observed for RE 
J03 1 7-853 - see abo ve - and [WR] cores of planetary neb- 
ulae PHamann 1997) ), and using the properties inferred for 
0.6 Mq remnant (L ~ L,„, ~ lO'^"^ gcm^ s"', R 



a 0.6 



10\m, ~ 0„ 



10 



-0.3, 



lO^'^^'-^gcm^s 
, we find f, 



m,sd 



10 Bg M_j yr, which is of the same order as the neutrino 
cooling timescale of ^ 10"* yr at the ignition line (for the whole 
range of remnants, 2 x 10^ < fmsd ^ 5 x lO'^yr). 

Accretion from the disk, loss of rotational support, and pos- 
sible cooling of the hot envelope could all compress and heat 
the remnant core. A detailed study of this is beyond the scope 
of this paper, but we can make first-order estimates of the ef- 



fects on our merger re mnants, and compare these with the 
more detailed analysis of Schwab et al. (2012 1 in one specific 
case. 

For our estimates, we make four assumptions: (i) spin- 
down and accretion are much faster than thermal processes, 
and do not lead to local dissipation (i.e., particles entropies are 
constant in time); (ii) all angular momentum is carried away 
to large distances; and (iii) corresponding matter ends up with 
zero total energy (i.e., is at large distances and has negligible 
kinetic and internal energy). From energy conservation, the 
last assumption implies that the remaining object will have 
the same total energy as our merger remnant (but a lower 
mass), the first that it will have the same entropy structure, and 
the second that it has no rotational support. To determine the 
properties, we first determine the entropy profile of the merger 
remnant, by averaging entropy over isopotential surfaces. We 
then use this entropy profile and an estimated central density 
to construct a spherically symmetric (non-spinning) hydro- 
static model, iterating on the central density until it has the 
correct total energy (inside of the zero-pressure surface). This 
automatically gives the mass contained in this object, which 
will be lower than our remnant mass, the remainder represent- 
ing material that, due to dissipation of rotational energy, has 
expanded out to large distances and therefore provides neg- 
ligible weight. To determine the evolution of hot spots, we 
order remnant particles by potential, and map them to their 
new positions in the final object, calculating new temperatures 
from the new densities, again assuming their entropy did not 
change (entropy is not constant over isopotential surfaces, so 
these temperatures are not strictly consistent with the hydro- 
static model). 



In Fig. 19 we show the results of our evolutionary estimate 
for our fiaucial 0.4 - 0.8 and 0.6 - 0.6 M© systems. For the 
former case, the core-envelope, originally 0.90 M©, accretes 
0.06 M© from the disk, the remaining 0.25 M© going to large 
distances. The central core is not significantly heated, while 
the lower-density hot envelope is, with the outer hot envelope 
along the rotational axis passing the ignition line. Since this 
material is almost non-degenerate, the resulting nuclear burn- 
ing will likely be stable (or be extinguished by expansion). 
Thus, not unexpected, the hot envelopes of dissimilar-mass 
mergers are not good candidates for a nuclear runaway. 

For the 0.6 - 0.6 M© system, the center of the final, spun- 
down object is at much higher density and temperature than 
the remnant, while much of the outer regions have become 
less dense and cool. The latter happens because similar-mass 
mergers have strong rotational support, and if this is removed 
their binding energy increases significantly. To compensate 
for this, a large amount of mass has to expand to large dis- 
tances, causing the core-envelope mass to decrease from 1.11 
Mq to 0.91 Mq. In the final object, the hottest point on the 
equatorial plane does not reach the ignition line, but the sig- 
nificantly hotter points above and below the equatorial plane 
do, at densities under which degeneracy pressure still domi- 
nates. Hence, if the hot spots indeed compress with the rest of 
the remnant, a thermonuclear runaway could be triggered. (Of 
course, a nuclear runaway would start as soon as the heating 
timescale becomes shorter than the compression timescale, 
which may happen closer to the ignition line.) 

In Fig. [18] we show the results of applying our estimates 
to all our merger remnants. One sees that all compress and 
heat, and almost every remnant whose accretor mass is above 
0.8 Mq will reach ignition somewhere on the equatorial plane, 
in many cases under degenerate conditions. We also chart 
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the evolution of the off-center hot spots in similar-mass merg- 
ers (square points), and while they are at lower density, they 
remain degenerate and are all pushed substantially further 
above the ignition line than their counterparts on the equato- 
rial plane. If these points can lead to nuclear runaways, then 
almost all similar-mass mergers with accretor mass above 0.5 
Mq could explode. 

The above suggests it is at least plausible that many of our 
mergers would eventually ignite in degenerate conditions, and 
that it thus is worthwhile to simulate their evolution in detail. 
Suitable simulations have recently been pioneered by! 



|et al. (2012) and Schwab et'lL] ( |20i2| l. |Shen et al. started with 
a one-dimensional simulation, where they porte d the remnant 
of a 0.6 - 0.9 Mq merger (from Dan et al.|201 1 1, and evolved 
it assuming a 7 = 5/3 polytropic equation or state and an 
a = 10"^ viscosity. They find the system spins down com- 
pletely due to outward angular momentum transport, and the 
rotationally-supported thick disk is transformed into a ten- 
uous, thermally-supported envelope that hardly affects the 
core. Over l onger, thermal evolu tion timescales (simulated 
using MESA, [Pax ton et al.|2011| l, this tenuous hot envelope 
cools, compresses the core, and lights off-center convective 
carbon burning, eventually turning the remnant into an ONe 
WD (that may end its l ife in an accretion induced collapse). 

[Schwab et al.| ( pOT2] l went a step further, porting the same 
0.6 - 0.9 Mq simulation, as well as seven other systems, into 
two-dimensional ZEUS-MP2 simulations ( [Hayes et al.|2006| l, 
using the Helmholtz equation of state and an a = 3 x 10"^ 
viscosity. They confirm the one-dimensional results, find- 
ing complete spin-down and transformation of the rotation- 
ally supported disk into a tenuous, spherically symmetric, hot 
envelope. They find a 50% increase in the temperature of the 
hottest point, and a factor of 3 increase in the corresponding 
density. They also find entropy to roughly be constant in the 
remnant, except in the outer regions and at the very center, 
where dissipation of rotational energy leads to heating. 

It is encouraging that the results of the above detailed sim- 
ulations are similar to what we find using our first-order es- 
timates. For our 0.6 - 0.9 Mq remnant, our estimate give in- 
creases for the hottest equatorial point of a factor of 2.5 in 
density a nd 1.6 in te rnperature. reasonably close to what is 
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found by Schwab et al!''! Thus, our simplifying assumptions 



appear to be appropriate at least for dissimilar-mass mergers, 
where most of the rotational dissipation will be in the disk and 
envelope, and the structure of the remainder is roughly spher- 
ically symmetric (both in density and temperature). It is not 
clear our estimates would be equally good for similar-mass 
mergers, where rotational dissipation should occur throughout 
the star, heating the entire remnant, and where there are sub- 
stantial differences between the remnant's density and tem- 
perature structures. It will thus be particularly interesting to 
simulate the further evolution of those remnants in more de- 



tail. U nf ortunately, no s uch remnants were included by Shen 
let al.l and lSchwab et al.l 

6.2. Possible Explosions? 



From our estimates, it seems that, as suggested by vKCJlO 
many merger remnants will ignite carbon fusion. If a detona 
tion is triggered, the resulting explosion may well resemble 

" For the hottest point along the rotational axis, we find a factor of 4.8 
increase in density and 2.1 increase in temperature, suggesting our model 
does not depict as well the evolution of the outer hot envelope along the 
rotational axis. 



an SN la. Indeed, if the remnants spun down before igni- 
tion, their structures ar e sufficiently close to that of a cold WD 
that the calculations of Sim et al. (2010 1 should apply. From 
our estimates, for mergers that have a total mass between 1 .2 
and 1.4M0 (which should be the most common ones), the fi- 
nal objects have masses between ^0.9 and ^1.1 Mq which 
matches fairly nicely the range of ~ 1 to ~ 1 .2 Mq required 



to reproduce the observed range of SN la luminosities ( Sim 
et al. 2010). 

Of course, it is far from clear whether ignition leads to a 
detonation, since we do not currently understand how deto 



nation s are triggered ( |Seitenzah l et al. 2009; Woosley et al. 
201 1 and references therein). Generally, it should help that 



Ignition in our remnants is at much lower density (a few 
10^ gcm"^) than is the case for near-Chandrasekhar models 



(^ 10 g cm" ), because complete b urning lea ds to much larger 
relative ove rpressures (e.g. Mazurek et al. [1977 ; SeitenzaU] 
[et al.|20(j9] l. Also, if a deflagration is starteoTplausible mech- 
anisms to transition to a detonation all seem to require densi- 
ties around 10^ gcm~^, where the conductive flame speeds are 
slower and the separation betw een the various burning fronts 
increases (e.g., [Woosley et aL|2 009 2011). 

An interesting aspect of our results is that for all cases ig- 
nition likely happens off-center: in shells for dissimilar-mass 
mergers and in hot spots along the rotational axis for similar- 
mass ones. Previous one-dimensional simulations suggested 
off-center ignition would lead to a slow deflagration flame 



that tu rns the CO WD into a ONe WD (e.g., Saio & Nomoto 
1985 I. However, these calculations assumed a hot spot many 
pressure scale heights above the center. For ignition closer 



to the cente r, a deflagration plume is produced (e. g., [Aspden 



etal. 


2011 


etal. 


2011 



, which may transition to a detonation ( Seitenzahl 
and unbind the star. 
Given our findings, it seems likely that, if a detonation 
occurs, it will be triggered off-center. It would be interest- 
ing to simulate the resulting explosion, and see whether one 
could reproduce the observational evidence for asymmetries, 
which have been interpreted in terms of off-center ignition 
(though so far only in the context of near-Chandrasekhar mod- 
els; Ma eda et al..2010 a b ). 

Finally, while we simulated only mergers of CO WDs, we 
can extrapolate our results to more massive ONe WDs. For 
these, the temperatures would be at least as high as for our 
1 Mq accretors, and, after further viscous evolution, the merg- 
ers should become hot enough to ignite Ne burning. If this 
also leads to a detonation, the lower fusion energy released 
would likely lead to a less energetic explosion than expected 
for a CO WD merger, but it would produce far more ^^Ni and 
have a very large mass. Plausibly, it would resemble an SN 
la like SN 2009dc, which had unusually low ejecta veloci- 
ties, produc ed ^ 1 .8Mq of ^^Ni and had a total ejecta mass of 
--2.8M0 ( |Taubenberger et al.|201 1) . 

7. CONCLUSION 

We have performed a large, detailed parameter-space study 
of CO WD mergers, extracting pertinent properties and pro- 
files for each remnant, and studying how these vary across pa- 
rameter space. For a merger involving dissimilar-mass WDs, 
with low qp = Pc,d/Pc,a> the outcome is a cold, slowly rotat- 
ing, degeneracy-supported remnant core, which is essentially 
unaffected by the merger, surrounded by a hot, roughly spher- 
ical envelope and, further out, by a sub-Keplerian disk. For 
a similar-mass merger, with high qp, an ellipsoidal core is 
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produced along with a small disk, and the entire remnant is 
hot and partly supported by rotation. The transition between 
these two regimes is smooth, but occurs roughly at qp ~ 0.6, 
or equivalently a mass difference AM = M^-Md ~ 0.1. We 
found that for a fixed qp, merger remnant curves are roughly 
homologous. We also presented trends for a number of merger 
remnant properties, providing linear scaling relations and best 
fits for most of them, hoping these can guide theoretical un- 
derstanding and help analytical estimates. 

We made first-order estimates of the post-merger viscous 
evolution and spindown, and found that it is plausible that a 
large fraction of the mergers simulated will even tually e xpe- 
rience a nuclear runaway, as was suggested by jVKCJlO and 
thus possibly end as thermonuclear supernovae. Further, de- 
tailed, simulations of this evolution across the whole param- 
eter sp ace, using techniqu es similar to those of Shen et aL] 
(12012) and Schwab et al.'('2012 ), would be required to con- 
firm this. If the evolution of these remnants results in a deto- 
nation, a detailed comparison of the resulting light curve with 
observations must be carried out. 

Our work represents one of the most detailed parameter 
studies of WD mergers to date. It would benefit, however, 
from resolution of a number of topics. First, for better preci- 
sion, it will be necessary to use better initial conditions. For 
synchronized systems, it is already known th is has nontrivial 
effect on the outcome ( |Dan et al.|2011|[2012[ ), and our results 
suggest it is important also for non-synchronized systems. 
Unfortunately, for the non-synchronized case, it is not trivial 
to implement the initial conditions, but better approximations 
are possible. Second, it would be useful to try to compare 
with merger simulations done with a grid code, which should 
have become more straightforward now that good moving 
mesh codes have b ecome available (Springel 2010; Duffell & 
|MacFadyen|201 1[ ). More generally, whether or not WDs are 
synchronized before the merger remains unknown, yet clearly 
affects the resulting merger. Hopefully, this can be resolved 
empirically, by measuring the spin frequency for WDs in the 
very short-period binarie s that have recently been discovered 
(e.g., [Brown etal.|2011 1. 
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9. APPENDIX 

We present CO WD binary pre-merger properties in Tab. [T] 
andp] and post-merger remnant properties in Tab.[3]-|7l 

Masses and distances are in Mq and Rq (Mq = 1.9891 x 
lO"*^ g, Rq = 6.95508 x 10'° cm), and all other quantities are 
in CGS. Numerical subscripts indicate powers of 10. 

Subscript "d" indicates donor, and "a" accretor, while 
sub/superscript "disk" and "ce" stand for the disk and core- 
envelope, respectively. Superscript "z" indicates a value along 
the rotational axis. Subscript "max" stands for the maxi- 



mum value of a property, either along the equatorial plane 
if standalone or along the rotational axis if paired with super- 
script "z". Superscript "cv" indicates an estimate of a prop- 
erty that takes into account post-merger convection along the 
equatorial plane (convection does not occur along the rota- 
tional axis). Subscript "rot" stands for rotational, "th" ther- 
mal, "deg" degeneracy, and "pot" gravitational potential en- 
ergy. 

The disk inner radius (Ddir (Tab. |7]i is determined by flatten- 
ing the system along the z-axis, and finding the iB at which 
the centripetal force counteracts half the gravitational force 
on the particles. Bound particles with (53 larger and smaller 
than this value are part of the disk and core-envelope, respec- 
tively (from this differentiation, disk and core-envelope bulk 
properties are determined). 

All bulk energies (E), moments of inertia (/) and angular 
momenta (L) are determined by summing up individual parti- 
cle values. Scaleheights are determined by fitting a Gaussian 
to the density profile of particles within a cylindrical shell of 
radius 03, except for h^^ which uses particles in the equato- 
rial plane. Enclosed masses (Menc) are determined spheri- 
cally (e.g. Menc(71nax) is the mass enclosed by a sphere of 
r = Ci!(7max))- In particular, energy enclosed masses Mi.^^(fE) 
are the spherical masses enclosing fraction / of energy E of 
the remnant (e.g. Menc(2£^th) is the spherical mass enclosing 
half of the total remnant thermal energy). 

In Tab. [3j Pc is determined by averaging values of particles 
within r < -j^/iz of the center This accounts for minor dis- 
crepancies between and Pmax. which is determined along 
the equatorial plane. Mgj is the total mass of particles that are 
unbound, which is determined by finding if their total energy 
£ = £po, + £K >0. 

In Tab. [4j central core mass Mcc represents the spherical 
mass enclosed within a radius r = h^, while Mcc.d represents 
the do nor mass enclosed; M^^ d/(-^cc~-^cc d) = (/d//a)cc (Sec. 
[3TT}. 

In Tab. |5J is determined in the same manner as pc, and 
therefore has discrepancies with the equatorially-determined 
Tmax when maximum temperature is at the center of the rem- 
nant. 

In Tab. [6] for some similar-mass mergers p(7',nax) will be 
larger than p(7max) {^{Tms.x) will also change); this is because 
after six orbital periods these mergers will have their maxi- 
mum temperature, but not necessarily their maximum density, 
exactly at the center of the remnant. We assume that the re- 
gion undergoing convection will evolve until it becomes isen- 
tropic; as a result maximum temperature after convection will 
shift to the point of maximum density. 

In Tab. |7] n3md is the cylindrical radius enclosing half the 
mass of the disk. 

The "0.625 - 0.65 255k(2)" simulation is a re-running of the 
parameter space 0.625 - 0.65 Mq run, to determine the extent 
to which random noise affects remnant properties (Sec. |4.7| i. 
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Table 1 

Table of WD binary initial conditions. 



Sim 


Mj 


Rd 


(Pc.d)b 


(^rot.t!)48 


(^th.d)48 


(^deg.d)48 


(^pot.ti)48 


(/d)51 




0.4 - 0.4 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.4 - 0.5 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.5 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.4 - 0.55 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.55 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55 - 0.55 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.4 - 0.6 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.6 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55 - 0.6 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.6 - 0.6 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.4 - 0.65 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.65 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55 - 0.65 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.575 - 0.65 


0.5780 


0.0116 


3.20 


0.00 


0.15 


49.27 


-91.30 


0.16 


0.0042 


0.6 - 0.65 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.625 - 0.65 


0.6285 


0.0110 


4.20 


0.00 


0.17 


62.67 


-114.37 


0.16 


0.0039 


0.64 - 0.65 


0.6399 


0.0108 


4.47 


0.00 


0.17 


66.05 


-120.10 


0.16 


0.0038 


0.65 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.4 - 0.7 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.7 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55 - 0.7 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.6 - 0.7 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.65 - 0.7 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.7 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.4 - 0.8 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.8 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55 - 0.8 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.6 - 0.8 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.65 - 0.8 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.7 - 0.8 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.8-0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.4 - 0.9 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5 - 0.9 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55 - 0.9 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.6 - 0.9 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.65 - 0.9 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.7 - 0.9 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.8 - 0.9 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.9 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.4 - 1.0 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.5-1.0 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.55-1.0 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.6-1.0 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.65-1.0 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.7- 1.0 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.8- 1.0 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.9-1.0 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


1.0-1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


Varying Composition 


0.4 - 0.4 He 


0.4023 


0.0145 


1.03 


0.00 


0.20 


17.90 


-34.67 


0.18 


0.0055 


0.4 - 0.8 He 


0.4023 


0.0145 


1.03 


0.00 


0.20 


17.90 


-34.67 


0.18 


0.0055 


0.4 - 0.4 Mg 


0.4023 


0.0134 


1.21 


0.00 


0.07 


18.57 


-36.77 


0.16 


0.0052 


0.4 - 0.8 Mg 


0.4023 


0.0134 


1.21 


0.00 


0.07 


18.57 


-36.77 


0.16 


0.0052 


Varying ao 


0.6 - 0.6 0.9a 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.6-0.6 l.lt! 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.6 - 0.8 0.9a 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.6-0.8 1.1a 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


Synchronization 


0.6 - 0.6 synch 


0.6035 


0.0114 


3.67 


1.42 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.6 - 0.8 synch 


0.6035 


0.0114 


3.67 


1.35 


0.16 


55.72 


-102.47 


0.16 


0.0040 


Varying Viscosity 


0.6 - 0.8 lowvisc 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.6 - 0.8 highvisc 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 








Varyin 


5 Numerical Resolution 










0.625 - 0.65 64k 


0.6282 


0.0109 


4.15 


0.00 


0.19 


62.31 


-113.81 


0.16 


0.0039 


0.625 - 0.65 128k 


0.6284 


0.0109 


4.18 


0.00 


0.17 


62.57 


-114.19 


0.16 


0.0039 


0.625 - 0.65 255k(2) 


0.6285 


0.0110 


4.20 


0.00 


0.17 


62.67 


-114.37 


0.16 


0.0039 


0.625 - 0.65 510k 


0.6284 


0.0111 


4.22 


0.00 


0.17 


62.70 


-114.40 


0.16 


0.0039 



Table 2 

Table of WD binary initial conditions, continued. 



Sim 






(Pc,a)6 


(^rot,a)48 


(^th,a)48 


(^deg,a)48 


(^pot.a)48 


(4)51 


Aa 


«o 


^"0* 


0.4 - 0.4 


0.4023 


0.0138 


1.14 


0.00 


0.09 


18.37 


-35.92 


0.17 


0.0053 


0.0364 


77.57 


0.4 ^ 0.5 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.0384 


79.05 


0.5 - 0.5 


0.5026 


0.0124 


2.08 


0.00 


0.12 


33.32 


-63.23 


0.17 


0.0046 


0.0328 


59.14 


0.4 - 0.55 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.0393 


79.61 


0.5 - 0.55 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.0335 


59.65 


0.55 - 0.55 


0.5533 


0.0118 


2.80 


0.00 


0.14 


43.53 


-81.29 


0.16 


0.0043 


0.0312 


52.42 


0.4 - 0.6 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0401 


80.07 


0.5 - 0.6 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0342 


60.08 


0.55 - 0.6 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0318 


52.83 


0.6 - 0.6 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0300 


47.29 


0.4 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0409 


80.46 


0.5 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0348 


60.45 


0.55 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0324 


53.19 


0.575 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0316 


50.57 


0.6 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0305 


47.63 


0.625 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0293 


44.35 


0.64 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0288 


42.90 


0.65 - 0.65 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0283 


41.67 


0.4 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.0416 


80.79 


0.5 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.0354 


60.77 


0.55 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.0330 


53.49 


0.6 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.03 1 1 


47.93 


0.65 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.0288 


41.95 


0.7 - 0.7 


0.7037 


0.0102 


6.28 


0.00 


0.20 


88.12 


-156.74 


0.15 


0.0035 


0.0268 


37.06 


0.4 - 0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0430 


81.31 


0.5 - 0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0366 


61.29 


0.55-0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0341 


54.00 


0.6-0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0321 


48.43 


0.65 - 0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0297 


42.42 


0.7 - 0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0277 


37.50 


0.8 - 0.8 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0242 


29.77 


0.4 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0443 


81.70 


0.5 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0377 


61.68 


0.55 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0351 


54.39 


0.6 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0330 


48.81 


0.65 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0306 


42.78 


0.7 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0284 


37.84 


0.8 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0249 


30.08 


0.9 - 0.9 


0.9049 


0.0082 


18.92 


0.00 


0.34 


203.49 


-335.60 


0.12 


0.0026 


0.0216 


23.66 


0.4-1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0456 


81.98 


0.5 - 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0387 


61.99 


0.55 - 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0360 


54.70 


0.6- 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0339 


49.12 


0.65 - 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0313 


43.08 


0.7 - 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0292 


38.13 


0.8 - 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0255 


30.34 


0.9- 1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0222 


23.88 


1.0-1.0 


1.0056 


0.0072 


34.93 


0.00 


0.45 


307.02 


-483.86 


0.09 


0.0022 


0.0191 


18.56 



Varying Composition 



0.4 - 0.4 He 


0.4023 


0.0145 


1.03 


0.00 


0.20 


17.90 


-34.67 


0.18 


0.0055 


0.0383 


83.44 


0.4 - 0.8 He 


0.8046 


0.0095 


9.82 


0.00 


0.46 


131.12 


-224.70 


0.14 


0.0031 


0.0452 


87.46 


0.4 - 0.4 Mg 


0.4023 


0.0134 


1.21 


0.00 


0.07 


18.57 


-36.77 


0.16 


0.0052 


0.0354 


74.35 


0.4 - 0.8 Mg 


0.8046 


0.0090 


11.65 


0.00 


0.23 


138.06 


-237.43 


0.13 


0.0030 


0.0418 


77.93 


Varying ao 


0.6 - 0.6 0.9a 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0270 


40.38 


0.6-0.6 1.1a 


0.6035 


0.0114 


3.67 


0.00 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0330 


54.56 


0.6 - 0.8 0.9a 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0289 


41.35 


0.6-0.8 1.1a 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0353 


55.87 


Synchronization 


0.6 - 0.6 synch 


0.6035 


0.0114 


3.67 


1.42 


0.16 


55.72 


-102.47 


0.16 


0.0040 


0.0300 


47.29 


0.6 - 0.8 synch 


0.8046 


0.0092 


10.84 


1.12 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0321 


48.43 


Varying Viscosity 


0.6 - 0.8 lowvisc 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0321 


48.43 


0.6 - 0.8 highvisc 


0.8046 


0.0092 


10.84 


0.00 


0.27 


135.13 


-231.87 


0.13 


0.0030 


0.0321 


48.43 



Varying Numerical Resolution 



0.625 - 0.65 64k 


0.6535 


0.0105 


4.72 


0.00 


0.20 


70.04 


-126.85 


0.16 


0.0038 


0.0292 


43.99 


0.625 - 0.65 128k 


0.6538 


0.0107 


4.78 


0.00 


0.19 


70.29 


-127.26 


0.16 


0.0038 


0.0289 


43.39 


0.625 - 0.65 255k(2) 


0.6536 


0.0107 


4.82 


0.00 


0.18 


70.40 


-127.39 


0.16 


0.0038 


0.0293 


44.35 


0.625 - 0.65 510k 


0.6537 


0.0108 


4.84 


0.00 


0.18 


70.50 


-127.55 


0.16 


0.0038 


0.0296 


45.08 



28 



Tables 

Table of merger remnant properties. 



Sim 


Time 


(Pc)6 


(Pmax)6 




ha 


Mce 






^disk,d 


Menc(jM<i) 


Me„c(iMd) 


Me„c(|Md) 


0.4 - 0.4 


466.3 


0.69 


0.67 


0.0044 


0.0092 


0.733 


0.000 


0.071 


0.036 


0.200 


0.403 


0.605 


0.4 - 0.5 


474.4 


2.13 


2.05 


0.0042 


0.0052 


0.700 


0.000 


0.205 


0.168 


0.512 


0.655 


0.790 


0.5 - 0.5 


354.3 


1.13 


1.13 


0.0039 


0.0085 


0.923 


0.001 


0.081 


0.044 


0.251 


0.506 


0.762 


0.4-0.55 


478.5 


2.87 


2.76 


0.0040 


0.0046 


0.726 


0.001 


0.229 


0.192 


0.572 


0.714 


0.845 


0.5 - 0.55 


358.4 


2.37 


2.29 


0.0038 


0.0059 


0.879 


0.002 


0.175 


0.144 


0.537 


0.733 


0.904 


0.55 - 0.55 


313.6 


1.44 


1.44 


0.0036 


0.0081 


1.018 


0.001 


0.087 


0.043 


0.280 


0.555 


0.832 


0.4 - 0.6 


480.5 


3.94 


3.78 


0.0038 


0.0041 


0.719 


0.001 


0.286 


0.251 


0.659 


0.777 


0.898 


0.5 - 0.6 


360.4 


3.59 


3.45 


0.0037 


0.0047 


0.876 


0.001 


0.229 


0.185 


0.608 


0.787 


0.960 


0.55-0.6 


317.6 


3.23 


3.11 


0.0036 


0.0054 


0.959 


0.002 


0.196 


0.160 


0.582 


0.801 


0.989 


0.6 - 0.6 


283.0 


1.77 


1.82 


0.0034 


0.0077 


1.108 


0.001 


0.098 


0.046 


0.303 


0.606 


0.906 


0.4 - 0.65 


482.6 


5.13 


4.91 


0.0036 


0.0038 


0.755 


0.001 


0.299 


0.272 


0.720 


0.835 


0.948 


0.5 - 0.65 


362.4 


4.82 


4.63 


0.0035 


0.0042 


0.891 


0.001 


0.264 


0.220 


0.674 


0.850 


1.019 


0.55 - 0.65 


319.7 


4.28 


4.11 


0.0034 


0.0048 


0.963 


0.001 


0.243 


0.197 


0.644 


0.851 


1.046 


0.575 - 0.65 


303.4 


3.98 


3.81 


0.0034 


0.0051 


1.002 


0.001 


0.229 


0.186 


0.625 


0.858 


1.058 


0.6 - 0.65 


285.1 


3.91 


3.75 


0.0034 


0.0053 


1.042 


0.001 


0.214 


0.170 


0.601 


0.856 


1.069 


0.625 - 0.65 


266.7 


3.82 


3.68 


0.0033 


0.0054 


1.072 


0.002 


0.208 


0.172 


0.571 


0.867 


1.093 


0.64 - 0.65 


256.6 


2.77 


2.67 


0.0033 


0.0066 


1.133 


0.001 


0.159 


0.131 


0.457 


0.816 


1.092 


0.65 - 0.65 


250.4 


2.34 


2.38 


0.0032 


0.0072 


1.190 


0.001 


0.117 


0.069 


0.343 


0.682 


1.013 


0.4 - 0.7 


484.6 


6.53 


6.26 


0.0034 


0.0036 


0.812 


0.001 


0.294 


0.273 


0.778 


0.892 


1.000 


0.5 - 0.7 


364.5 


6.46 


6.22 


0.0033 


0.0037 


0.889 


0.001 


0.316 


0.265 


0.746 


0.909 


1.071 


0.55 - 0.7 


321.7 


6.23 


5.98 


0.0033 


0.0039 


0.978 


0.001 


0.277 


0.231 


0.725 


0.917 


1.105 


0.6 - 0.7 


287.1 


5.69 


5.48 


0.0032 


0.0044 


1.047 


0.001 


0.259 


0.213 


0.696 


0.923 


1.134 


0.65 - 0.7 


252.5 


5.13 


4.92 


0.0031 


0.0049 


1.129 


0.002 


0.227 


0.182 


0.671 


0.932 


1.155 


0.7 - 0.7 


221.9 


2.90 


2.99 


0.0031 


0.0068 


1.285 


0.001 


0.122 


0.065 


0.359 


0.715 


1.070 


0.4 - 0.8 


488.7 


10.93 


10.46 


0.0030 


0.0031 


0.900 


0.001 


0.306 


0.296 


0.895 


1.000 


1.103 


0.5 - 0.8 


368.5 


11.56 


11.09 


0.0029 


0.0031 


0.927 


0.001 


0.379 


0.346 


0.891 


1.033 


1.173 


0.55 - 0.8 


323.7 


11.99 


11.51 


0.0029 


0.0031 


0.945 


0.002 


0.411 


0.368 


0.889 


1.048 


1.209 


0.6 - 0.8 


291.2 


11.67 


11.21 


0.0028 


0.0031 


1.011 


0.001 


0.395 


0.333 


0.865 


1.051 


1.243 


0.65 - 0.8 


254.5 


11.09 


10.68 


0.0028 


0.0033 


1.117 


0.002 


0.339 


0.280 


0.833 


1.055 


1.276 


0.7 - 0.8 


224.0 


9.82 


9.44 


0.0028 


0.0038 


1.208 


0.002 


0.298 


0.239 


0.806 


1.054 


1.301 


0.8 - 0.8 


179.2 


4.31 


4.55 


0.0027 


0.0063 


1.489 


0.001 


0.119 


0.059 


0.410 


0.809 


1.208 


0.4 - 0.9 


490.7 


18.45 


17.64 


0.0026 


0.0027 


0.988 


0.001 


0.319 


0.313 


1.001 


1.103 


1.205 


0.5 - 0.9 


370.6 


19.59 


18.81 


0.0025 


0.0026 


1.021 


0.001 


0.385 


0.370 


1.014 


1.147 


1.277 


0.55 - 0.9 


325.8 


20.15 


19.30 


0.0025 


0.0026 


1.027 


0.001 


0.430 


0.403 


1.014 


1.166 


1.312 


0.6 - 0.9 


293.2 


21.07 


20.19 


0.0025 


0.0026 


1.060 


0.002 


0.446 


0.408 


1.010 


1.179 


1.346 


0.65 - 0.9 


256.6 


21.66 


20.75 


0.0024 


0.0026 


1.060 


0.002 


0.497 


0.444 


1.005 


1.192 


1.382 


0.7 - 0.9 


228.0 


20.59 


19.74 


0.0024 


0.0027 


1.164 


0.002 


0.442 


0.372 


0.970 


1.191 


1.417 


0.8 - 0.9 


181.2 


18.61 


17.86 


0.0024 


0.0031 


1.333 


0.003 


0.373 


0.305 


0.932 


1.204 


1.478 


0.9 - 0.9 


142.5 


8.45 


8.16 


0.0024 


0.0052 


1.614 


0.002 


0.194 


0.138 


0.525 


1.029 


1.473 


0.4 - 1.0 


492.7 


32.25 


30.90 


0.0022 


0.0023 


1.080 


0.002 


0.326 


0.323 


1.103 


1.205 


1.306 


0.5 - 1.0 


372.6 


34.15 


32.67 


0.0022 


0.0022 


1.100 


0.001 


0.407 


0.399 


1.124 


1.252 


1.380 


0.55 - 1.0 


327.8 


35.34 


33.89 


0.0021 


0.0022 


1.116 


0.001 


0.441 


0.427 


1.131 


1.274 


1.416 


0.6-1.0 


295.2 


36.65 


35.21 


0.0021 


0.0022 


1.132 


0.002 


0.475 


0.455 


1.136 


1.296 


1.452 


0.65-1.0 


258.6 


38.54 


37.04 


0.0021 


0.0022 


1.173 


0.003 


0.483 


0.451 


1.131 


1.311 


1.485 


0.7- 1.0 


228.0 


40.51 


38.87 


0.0020 


0.0022 


1.177 


0.003 


0.529 


0.485 


1.130 


1.325 


1.521 


0.8- 1.0 


181.2 


40.89 


39.19 


0.0020 


0.0022 


1.315 


0.003 


0.492 


0.423 


1.096 


1.342 


1.590 


0.9-1.0 


142.5 


39.41 


37.75 


0.0019 


0.0023 


1.457 


0.004 


0.449 


0.367 


1.055 


1.349 


1.648 


1.0-1.0 


112.0 


27.94 


26.81 


0.0019 


0.0031 


1.625 


0.007 


0.379 


0.071 


0.312 


0.643 


1.060 


Varying Composition 


0.4 - 0.4 He 


466.3 


0.56 


0.57 


0.0046 


0.0099 


0.736 


0.000 


0.068 


0.038 


0.203 


0.411 


0.614 


0.4 - 0.8 He 


488.7 


9.74 


9.34 


0.0031 


0.0032 


0.894 


0.000 


0.313 


0.304 


0.897 


1.001 


1.104 


0.4 - 0.4 Mg 


466.3 


0.75 


0.73 


0.0043 


0.0090 


0.737 


0.000 


0.068 


0.032 


0.199 


0.397 


0.599 


0.4 - 0.8 Mg 


488.7 


11.81 


11.33 


0.0029 


0.0030 


0.906 


0.001 


0.300 


0.290 


0.893 


1.000 


1.103 


Varying oq 


0.6 -0.6 0.9a 


283.0 


2.03 


2.01 


0.0034 


0.0073 


1.109 


0.000 


0.098 


0.047 


0.300 


0.604 


0.905 


0.6-0.6 1.1a 


724.9 


2.23 


2.13 


0.0035 


0.0072 


1.045 


0.001 


0.162 


0.127 


0.404 


0.728 


1.007 


0.6 - 0.8 0.9a 


230.1 


11.55 


11.08 


0.0028 


0.0032 


1.059 


0.002 


0.347 


0.288 


0.843 


1.042 


1.235 


0.6-0.8 l.lo 


861.3 


11.87 


11.36 


0.0028 


0.0031 


0.984 


0.002 


0.422 


0.374 


0.892 


1.067 


1.243 


Synchronization 


0.6 - 0.6 synch 


407.2 


2.96 


2.84 


0.0034 


0.0060 


1.002 


0.001 


0.204 


0.070 


0.268 


0.544 


0.831 


0.6 - 0.8 synch 


313.6 


11.58 


11.11 


0.0028 


0.0032 


0.984 


0.002 


0.422 


0.384 


0.907 


1.077 


1.247 


Varying Viscosity 


0.6 - 0.8 lowvisc 


197.5 


11.92 


11.39 


0.0028 


0.0031 


0.968 


0.000 


0.440 


0.370 


0.852 


1.050 


1.236 


0.6 - 0.8 highvisc 


242.3 


11.30 


10.88 


0.0028 


0.0032 


1.018 


0.001 


0.389 


0.328 


0.864 


1.052 


1.243 


Varying Numerical Resolution 


0.625 - 0.65 64k 


266.7 


3.93 


3.78 


0.0033 


0.0054 


1.067 


0.002 


0.212 


0.173 


0.551 


0.852 


1.092 


0.625 - 0.65 128k 


266.7 


3.54 


3.39 


0.0033 


0.0057 


1.091 


0.002 


0.190 


0.159 


0.533 


0.857 


1.095 


0.625 - 0.65 255k(2) 


266.7 


3.74 


3.59 


0.0033 


0.0056 


1.077 


0.002 


0.203 


0.167 


0.573 


0.869 


1 .090 


0.625 -0.65 510k 


266.7 


3.18 


3.05 


0.0033 


0.0060 


1.103 


0.002 


0.177 


0.138 


0.573 


0.867 


1.080 



Table 4 

Table of merger remnant properties, continued. 



Sim 






(^rot)48 


(£th)48 


(^deg)48 


(^pot)48 


(£rof)48 


(£ft*)48 




(ES)48 


(^uf)48 




0.4 


0.4 


0.038 


0.019 


14.58 


4.56 


25.94 


-88.53 


1.97 


0.42 


0.10 


12.61 


4.15 


25.84 


0.4 


0.5 


0.090 


0.000 


14.82 


7.26 


40.58 


-120.04 


6.31 


1.48 


0.65 


8.51 


5.78 


39.94 


0.5 


0.5 


0.043 


0.022 


25.94 


7.47 


46.12 


-153.92 


3.18 


0.66 


0.17 


22.77 


6.82 


45.95 


0.4- 


0.55 


0.102 


0.000 


15.58 


8.02 


50.35 


-139.92 


7.65 


1.93 


0.79 


7.93 


6.09 


49.56 


0.5 - 


0.55 


0.077 


0.003 


25.06 


10.63 


56.48 


-175.68 


7.09 


1.58 


0.63 


17.97 


9.05 


55.85 


0.55 


0.55 


0.046 


0.023 


33.41 


9.29 


59.42 


-196.19 


3.92 


0.79 


0.22 


29.48 


8.50 


59.20 


0.4 


0.6 


0.119 


0.000 


15.48 


8.62 


63.53 


-163.25 


10.61 


2.72 


1.38 


4.86 


5.91 


62.15 


0.5 


0.6 


0.100 


0.000 


25.21 


12.41 


68.81 


-199.98 


9.12 


2.27 


0.85 


16.09 


10.14 


67.96 


0.55 


-0.6 


0.085 


0.003 


31.57 


13.52 


72.89 


-222.47 


9.26 


2.14 


0.90 


22.31 


11.38 


72.00 


0.6 


0.6 


0.049 


0.024 


42.01 


11.30 


75.14 


-244.95 


5.04 


0.99 


0.31 


36.97 


10.31 


74.83 


0.4- 


0.65 


0.130 


0.000 


15.72 


9.58 


78.30 


-190.08 


11.19 


2.98 


1.23 


4.52 


6.60 


iim 


0.5- 


0.65 


0.115 


0.000 


26.24 


13.60 


82.98 


-227.54 


11.04 


3.02 


1.14 


15.20 


10.58 


81.83 


0.55 


0.65 


0.098 


0.001 


32.83 


15.38 


85.92 


-250.14 


11.82 


3.03 


1.24 


21.01 


12.35 


84.68 


0.575 


-0.65 


0.090 


0.004 


36.32 


16.03 


87.91 


-262.10 


11.12 


2.70 


1.04 


25.21 


13.33 


86.86 


0.6- 


0.65 


0.087 


0.005 


39.89 


16.45 


90.80 


-275.67 


10.32 


2.42 


0.96 


29.57 


14.03 


89.84 


0.625 


-0.65 


0.083 


0.016 


43.20 


17.02 


95.69 


-291.23 


11.78 


2.80 


1.21 


31.41 


14.22 


94.48 


0.64 


0.65 


0.063 


0.019 


48.76 


15.03 


92.84 


-295.23 


9.54 


1.83 


0.89 


39.22 


13.20 


91.95 


0.65 


0.65 


0.054 


0.026 


51.55 


14.34 


94.43 


-302.89 


7.12 


1.43 


0.54 


44.43 


12.91 


93.89 


0.4 


0.7 


0.139 


0.000 


16.37 


10.50 


95.10 


-220.55 


11.67 


3.06 


1.00 


4.69 


7.44 


94.09 


0.5 


0.7 


0.128 


0.000 


27.25 


14.65 


100.20 


-259.34 


13.46 


3.73 


1.50 


13.78 


10.91 


98.70 


0.55 


-0.7 


0.120 


0.000 


33.05 


17.10 


104.47 


-283.36 


12.76 


3.60 


1.23 


20.29 


13.50 


103.24 


0.6 


0.7 


0.107 


0.001 


40.42 


18.78 


108.42 


-309.13 


13.64 


3.56 


1.36 


26.78 


15.22 


107.06 


0.65 


-0.7 


0.093 


0.005 


48.92 


20.58 


1 14.00 


-339.51 


13.55 


3.22 


1.37 


35.37 


17.37 


112.63 


0.7 


0.7 


0.056 


0.027 


63.00 


17.36 


116.76 


-369.34 


8.16 


1.71 


0.64 


54.84 


15.65 


116.13 


0.4 


0.8 


0.156 


0.000 


17.50 


12.53 


140.67 


-298.14 


13.53 


3.47 


0.90 


3.97 


9.06 


139.77 


0.5 


0.8 


0.152 


0.000 


27.26 


17.07 


149.82 


-341.15 


19.70 


5.54 


2.28 


7.57 


11.53 


147.54 


0.55 


-0.8 


0.151 


0.000 


32.96 


19.43 


156.29 


-367.57 


23.43 


6.52 


3.46 


9.53 


12.91 


152.83 


0.6 


0.8 


0.144 


0.000 


41.48 


22.64 


158.16 


-394.06 


22.63 


6.18 


2.93 


18.84 


16.46 


155.23 


0.65 


-0.8 


0.135 


0.000 


49.52 


25.83 


163.48 


-426.06 


19.76 


5.65 


2.18 


29.76 


20.18 


161.29 


0.7 


0.8 


0.118 


0.001 


59.22 


29.01 


168.15 


-460.68 


20.03 


5.64 


2.22 


39.19 


23.37 


165.93 


0.8 


0.8 


0.060 


0.029 


91.69 


24.23 


173.98 


-533.98 


9.84 


2.13 


0.65 


81.85 


22.10 


173.32 


0.4 


0.9 


0.171 


0.000 


18.43 


15.42 


205.64 


-402.64 


14.53 


3.74 


0.82 


3.90 


11.69 


204.83 


0.5 


0.9 


0.169 


0.000 


28.95 


19.77 


216.90 


-449.00 


22.11 


5.77 


1.79 


6.84 


14.00 


215.12 


0.55 


-0.9 


0.168 


0.000 


35.58 


22.61 


222.87 


-476.80 


26.50 


7.07 


2.67 


9.08 


15.54 


220.19 


0.6 


0.9 


0.166 


0.000 


40.82 


26.24 


232.20 


-508.73 


29.47 


8.83 


4.11 


11.35 


17.42 


228.09 


0.65 


-0.9 


0.164 


0.000 


49.12 


30.14 


240.24 


-544.63 


34.66 


10.36 


5.58 


14.46 


19.78 


234.66 


0.7 


0.9 


0.156 


0.000 


60.85 


34.53 


241.46 


-579.54 


31.53 


9.54 


4.31 


29.32 


24.99 


237.16 


0.8 


0.9 


0.136 


0.000 


82.83 


42.72 


258.55 


-668.64 


30.10 


8.93 


3.74 


52.74 


33.79 


254.81 


0.9 


0.9 


0.076 


0.032 


125.32 


36.34 


257.45 


-756.81 


20.12 


4.42 


2.28 


105.19 


31.92 


255.17 


0.4 


1.0 


0.185 


0.000 


19.09 


19.37 


300.71 


-546.26 


15.45 


4.72 


0.74 


3.65 


14.65 


299.97 


0.5 


1.0 


0.183 


0.000 


31.08 


23.85 


313.10 


-596.36 


24.90 


6.57 


1.79 


6.18 


17.28 


311.31 


0.55 


- 1.0 


0.182 


0.000 


37.66 


26.93 


321.49 


-627.22 


29.61 


8.06 


2.47 


8.05 


18.87 


319.02 


0.6 


1.0 


0.181 


0.000 


44.78 


29.81 


330.80 


-660.58 


34.50 


9.55 


3.47 


10.28 


20.26 


327.34 


0.65 


- 1.0 


0.179 


0.000 


50.67 


35.72 


344.54 


-702.78 


36.41 


11.85 


4.60 


14.26 


23.87 


339.94 


0.7 


1.0 


0.178 


0.000 


58.86 


40.62 


358.64 


-748.26 


41.51 


14.00 


6.42 


17.35 


26.62 


352.22 


0.8 


1.0 


0.169 


0.000 


82.52 


52.77 


376.89 


-844.96 


42.44 


14.76 


5.96 


40.08 


38.02 


370.93 


0.9 


1.0 


0.157 


0.000 


111.99 


63.23 


402.93 


-964.88 


45.25 


14.73 


6.49 


66.74 


48.50 


396.44 


1.0 


1.0 


0.119 


0.096 


156.75 


68.47 


423.55 


-1103.35 


49.67 


13.45 


8.16 


107.08 


55.02 


415.39 



Varying Composition 



0.4 - 0.4 He 


0.035 


0.018 


14.10 


4.30 


25.00 -84.66 


1.78 


0.37 


0.10 


12.32 


3.94 


24.90 


0.4 - 0.8 He 


0.157 


0.000 


17.13 


11.38 


134.62 -285.74 


13.80 


3.15 


1.10 


3.34 


8.23 


133.52 


0.4 - 0.4 Mg 


0.039 


0.020 


15.02 


4.68 


26.25 -90.97 


1.92 


0.41 


0.08 


13.09 


4.27 


26.17 


0.4 - 0.8 Mg 


0.155 


0.000 


17.79 


13.23 


144.38 -306.53 


13.39 


3.63 


0.79 


4.41 


9.59 


143.58 


Varying oq 


0.6 - 0.6 0.9a 


0.054 


0.027 


41.21 


13.46 


77.86 -251.57 


5.21 


1.18 


0.39 


35.99 


12.28 


77.47 


0.6-0.6 1.1a 


0.058 


0.020 


40.14 


10.78 


74.57 -240.74 


8.12 


1.33 


0.69 


32.02 


9.45 


73.88 


0.6 - 0.8 0.9a 


0.141 


0.000 


42.10 


24.59 


159.82 -401.75 


19.49 


5.85 


2.30 


22.61 


18.73 


157.52 


0.6-0.8 1.1a 


0.148 


0.000 


38.86 


21.06 


159.09 -387.31 


26.50 


7.04 


4.04 


12.36 


14.02 


155.05 


Synchronization 


0.6 - 0.6 synch 


0.073 


0.041 


39.95 


8.95 


77.63 -240.10 


9.69 


1.85 


0.94 


30.27 


7.10 


76.69 


0.6 - 0.8 synch 


0.143 


0.000 


42.06 


17.23 


158.47 -386.83 


27.99 


6.65 


4.42 


14.06 


10.58 


154.06 


Varying Viscosity 


0.6 - 0.8 lowvisc 


0.148 


0.000 


42.97 


21.36 


157.51 -394.55 


30.40 


6.37 


4.87 


12.57 


14.98 


152.64 


0.6 - 0.8 highvisc 


0.141 


0.000 


41.83 


22.75 


155.88 -392.60 


22.83 


5.92 


3.07 


19.00 


16.83 


152.81 



Varying Numerical Resolution 



0.625 - 0.65 64k 


0.084 


0.013 


43.14 


17.25 


95.31 


-290.86 


11.45 


2.86 


1.24 


31.69 


14.40 


94.08 


0.625 - 0.65 128k 


0.078 


0.015 


44.63 


16.60 


93.97 


-290.86 


10.90 


2.51 


1.11 


33.73 


14.09 


92.86 


0.625 - 0.65 255k(2) 


0.081 


0.013 


43.52 


16.77 


95.07 


-290.84 


11.59 


2.67 


1.17 


31.93 


14.10 


93.89 


0.625 - 0.65 510k 


0.072 


0.014 


44.92 


16.17 


92.73 


-288.83 


10.35 


2.14 


0.94 


34.58 


14.03 


91.79 
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Tables 

Table of merger remnant properties, continued. 



Sim 


(7c)8 (7max)8 


ra(rmax) {p{Tmsa))6 {S{Tmm)h Menc{Tm!a) {T^^h 


(p(7Sax))6 M 




) Menc(j£th) M 


encC^fith) M 


enc( 4^th) 


0.4 - 0.4 


1.29 


1.48 


0.0003 


0.66 


1.19 


0.000 


1.95 


0.23 


0.116 


0.199 


0.400 


0.608 


0.4 - 0.5 


0.55 


1.74 


0.0092 


0.43 


1.36 


0.424 


2.01 


0.08 


0.501 


0.403 


0.542 


0.681 


0.5 - 0.5 


1.86 


2.03 


0.0005 


1.09 


1.20 


0.000 


2.71 


0.39 


0.125 


0.243 


0.497 


0.761 


0.4 - 0.55 


0.37 


1.97 


0.0090 


0.44 


1.41 


0.484 


2.35 


0.11 


0.531 


0.471 


0.597 


0.723 


0.5 - 0.55 


1.03 


1.58 


0.0101 


0.53 


1.26 


0.577 


2.58 


0.24 


0.418 


0.366 


0.602 


0.807 


0.55 - 0.55 


1.97 


2.40 


0.0004 


1.37 


1.22 


0.000 


3.16 


0.57 


0.110 


0.264 


0.544 


0.835 


0.4 - 0.6 


0.39 


2.47 


0.0090 


0.38 


1.57 


0.557 


2.70 


0.07 


0.628 


0.546 


0.647 


0.751 


0.5 - 0.6 


0.82 


2.16 


0.0083 


0.76 


1.31 


0.517 


2.77 


0.16 


0.583 


0.463 


0.657 


0.838 


0.55 - 0.6 


0.96 


1.87 


0.0094 


0.69 


1.27 


0.624 


2.91 


0.31 


0.474 


0.411 


0.660 


0.884 


0.6 - 0.6 


2.71 


2.92 


0.0003 


1.70 


1.25 


0.000 


3.60 


0.81 


0.095 


0.283 


0.590 


0.912 


0.4 - 0.65 


0.28 


2.87 


0.0088 


0.37 


1.67 


0.619 


3.01 


0.06 


0.692 


0.600 


0.693 


0.786 


0.5 - 0.65 


0.83 


2.50 


0.0082 


0.75 


1.38 


0.582 


3.09 


0.10 


0.691 


0.538 


0.712 


0.879 


0.55 - 0.65 


1.04 


2.29 


0.0087 


0.81 


1.32 


0.635 


3.09 


0.16 


0.645 


0.472 


0.712 


0.922 


0.575 - 0.65 


1.28 


2.06 


0.0092 


0.77 


1.28 


0.689 


3.30 


0.51 


0.452 


0.432 


0.708 


0.942 


0.6 - 0.65 


1.47 


2.02 


0.0082 


1.11 


1.20 


0.594 


3.31 


0.45 


0.460 


0.421 


0.707 


0.962 


0.625 - 0.65 


2.11 


2.29 


0.0003 


3.67 


1.04 


0.000 


3.50 


0.63 


0.395 


0.415 


0.729 


0.997 


0.64 - 0.65 


2.34 


2.77 


0.0004 


2.65 


1.14 


0.000 


4.00 


0.72 


0.234 


0.343 


0.676 


1.002 


0.65 - 0.65 


3.11 


3.30 


0.0004 


2.24 


1.24 


0.000 


4.18 


0.91 


0.134 


0.312 


0.650 


0.998 


0.4 - 0.7 


0.85 


3.22 


0.0087 


0.36 


1.75 


0.677 


3.42 


0.07 


0.738 


0.650 


0.739 


0.824 


0.5 - 0.7 


0.55 


2.96 


0.0078 


0.79 


1.45 


0.625 


3.55 


0.16 


0.707 


0.619 


0.763 


0.914 


0.55 - 0.7 


1.14 


2.76 


0.0077 


1.01 


1.35 


0.623 


3.47 


0.22 


0.690 


0.572 


0.769 


0.958 


0.6 - 0.7 


1.43 


2.51 


0.0078 


1.16 


1.27 


0.645 


3.52 


0.29 


0.646 


0.524 


0.776 


1.001 


0.65 - 0.7 


1.77 


2.28 


0.0084 


1.17 


1.23 


0.715 


3.90 


0.47 


0.549 


0.464 


0.772 


1.042 


0.7 - 0.7 


3.75 


3.78 


0.0002 


2.79 


1.25 


0.000 


4.78 


0.95 


0.181 


0.334 


0.701 


1.076 


0.4 - 0.8 


1.00 


4.05 


0.0082 


0.35 


1.93 


0.789 


4.15 


0.30 


0.776 


0.755 


0.832 


0.898 


0.5 - 0.8 


0.54 


4.18 


0.0074 


0.72 


1.68 


0.760 


4.35 


0.11 


0.850 


0.743 


0.855 


0.969 


0.55 - 0.8 


0.35 


4.27 


0.0073 


0.83 


1.64 


0.762 


4.53 


0.15 


0.851 


0.740 


0.867 


1.006 


0.6 - 0.8 


0.50 


4.09 


0.0069 


1.19 


1.50 


0.729 


4.65 


0.30 


0.801 


0.720 


0.877 


1.053 


0.65 - 0.8 


1.21 


3.72 


0.0066 


1.68 


1.36 


0.702 


4.56 


0.31 


0.810 


0.671 


0.888 


1.106 


0.7 - 0.8 


1.67 


3.34 


0.0074 


1.60 


1.32 


0.800 


4.64 


0.99 


0.606 


0.595 


0.890 


1.150 


0.8-0.8 


4.82 


4.94 


0.0002 


4.13 


1.28 


0.000 


5.95 


1.97 


0.118 


0.369 


0.790 


1.226 


0.4 - 0.9 


1.56 


5.36 


0.0076 


0.35 


2.15 


0.895 


5.49 


0.30 


0.890 


0.853 


0.924 


0.981 


0.5 - 0.9 


1.50 


5.17 


0.0070 


0.70 


1.83 


0.879 


5.27 


0.12 


0.949 


0.851 


0.946 


1.033 


0.55 - 0.9 


1.13 


5.30 


0.0067 


0.94 


1.74 


0.867 


5.42 


0.17 


0.949 


0.845 


0.957 


1.067 


0.6 - 0.9 


1.33 


5.49 


0.0065 


1.17 


1.69 


0.861 


5.68 


0.29 


0.933 


0.843 


0.969 


1.110 


0.65 - 0.9 


0.86 


5.71 


0.0061 


1.57 


1.61 


0.838 


5.78 


0.35 


0.935 


0.834 


0.981 


1.148 


0.7 - 0.9 


0.60 


5.15 


0.0059 


2.23 


1.45 


0.798 


5.85 


0.47 


0.913 


0.803 


0.998 


1.207 


0.8 - 0.9 


2.18 


4.48 


0.0065 


2.45 


1.35 


0.908 


5.97 


1.23 


0.789 


0.710 


1.016 


1.297 


0.9 - 0.9 


4.42 


4.65 


0.0002 


8.10 


1.13 


0.000 


6.88 


2.04 


0.314 


0.467 


0.950 


1.422 


0.4 - 1.0 


1.49 


6.81 


0.0069 


0.39 


2.34 


1.000 


6.89 


0.34 


0.996 


0.917 


1.016 


1.067 


0.5 - 1.0 


2.32 


6.55 


0.0066 


0.65 


2.06 


0.996 


6.79 


0.56 


0.987 


0.947 


1.034 


1.107 


0.55 - 1.0 


2.19 


6.60 


0.0063 


0.92 


1.91 


0.987 


6.65 


0.65 


0.988 


0.947 


1.046 


1.135 


0.6-1.0 


1.94 


6.59 


0.0060 


1.26 


1.79 


0.974 


6.89 


0.25 


1.047 


0.948 


1.059 


1.168 


0.65-1.0 


1.89 


7.11 


0.0056 


1.77 


1.72 


0.957 


7.10 


0.31 


1.058 


0.944 


1.073 


1.214 


0.7-1.0 


1.35 


7.36 


0.0055 


2.03 


1.70 


0.960 


7.27 


0.53 


1.039 


0.940 


1.085 


1.256 


0.8-1.0 


1.14 


7.06 


0.0050 


3.76 


1.48 


0.900 


7.71 


1.49 


0.948 


0.905 


1.116 


1.349 


0.9 - 1.0 


2.34 


6.34 


0.0048 


5.42 


1.33 


0.890 


7.75 


1.75 


0.936 


0.835 


1.138 


1.437 


1.0-1.0 


6.38 


6.26 


0.0002 


26.80 


1.04 


0.000 


9.43 


7.41 


0.406 


0.678 


1.171 


1.556 


Varying Composition 


0.4 - 0.4 He 


0.76 


0.95 


0.0005 


0.54 


2.72 


0.000 


1.36 


0.20 


0.099 


0.189 


0.391 


0.602 


0.4 - 0.8 He 


0.25 


2.92 


0.0090 


0.22 


3.89 


0.800 


3.25 


0.03 


0.852 


0.754 


0.826 


0.891 


0.4 - 0.4 Mg 


1.58 


1.68 


0.0003 


0.72 


0.85 


0.000 


2.12 


0.24 


0.122 


0.202 


0.403 


0.608 


0.4 - 0.8 Mg 


1.52 


4.47 


0.0080 


0.39 


1.54 


0.790 


4.52 


0.32 


0.778 


0.756 


0.834 


0.902 


Varying oq 


0.6 -0.6 0.9a 


2.76 


2.91 


0.0004 


1.95 


1.22 


0.000 


3.90 


0.50 


0.226 


0.293 


0.598 


0.918 


0.6-0.6 l.lti 


1.33 


1.72 


0.0007 


2.12 


1.04 


0.001 


3.14 


0.57 


0.220 


0.341 


0.648 


0.933 


0.6 - 0.8 0.9a 


0.64 


4.12 


0.0067 


1.42 


1.45 


0.709 


4.39 


0.41 


0.776 


0.689 


0.871 


1.065 


0.6-0.8 1.1a 


0.97 


3.98 


0.0071 


1.01 


1.54 


0.757 


4.74 


0.21 


0.836 


0.734 


0.881 


1.035 


Synchronization 


0.6 - 0.6 synch 


0.89 


1.32 


0.0104 


0.65 


1.16 


0.720 


2.32 


0.18 


0.651 


0.416 


0.703 


0.973 


0.6 - 0.8 synch 


0.63 


3.50 


0.0074 


0.95 


1.48 


0.778 


4.26 


0.15 


0.846 


0.765 


0.902 


1.061 


Varying Viscosity 


0.6 - 0.8 lowvisc 


0.42 


4.22 


0.0068 


1.20 


1.52 


0.713 


4.12 


0.27 


0.816 


0.683 


0.844 


0.998 


0.6 - 0.8 highvisc 


0.36 


4.45 


0.0069 


1.18 


1.55 


0.723 


4.69 


0.38 


0.776 


0.723 


0.871 


1.044 


Varying Numerical Resolution 


0.625 - 0.65 64k 


1.27 


2.14 


0.0084 


1.11 


1.22 


0.628 


3.52 


0.56 


0.423 


0.434 


0.728 


0.997 


0.625 - 0.65 128k 


1.63 


2.08 


0.0005 


3.38 


1.02 


0.000 


3.44 


0.93 


0.272 


0.394 


0.706 


0.996 


0.625 - 0.65 255k(2) 


1.87 


2.15 


0.0004 


3.57 


1.02 


0.000 


3.56 


0.39 


0.476 


0.406 


0.721 


0.999 


0.625 -0.65 510k 


2.78 


2.79 


0.0003 


3.05 


1.12 


0.000 


3.89 


0.91 


0.238 


0.366 


0.702 


0.991 



Table 6 

Table of merger remnant properties, continued. 



Sim {T^^h aSiT^) (p(r^'ax))6 {S{TZx)h He Omax aJCHmax) M^nci^n,^) MencC^firot) Me„o(^£rot) Me„c(|£rot) 



0.4 - 0.4 


0.96 


0.0031 


0.67 


1.05 


0.35 


0.35 


0.0003 


0.000 


0.315 


0.515 


0.675 


0.4 - 0.5 


1.74 


0.0092 


0.43 


1.36 


0.16 


0.22 


0.0104 


0.489 


0.532 


0.668 


0.770 


0.5 - 0.5 


1.24 


0.0034 


1.13 


1.04 


0.46 


0.46 


0.0003 


0.000 


0.387 


0.635 


0.839 


0.4 - 0.55 


1.97 


0.0090 


0.44 


1.41 


0.15 


0.24 


0.0104 


0.555 


0.589 


0.713 


0.810 


0.5 - 0.55 


1.58 


0.0101 


0.53 


1.26 


0.34 


0.34 


0.0003 


0.000 


0.516 


0.742 


0.896 


0.55 - 0.55 


1.37 


0.0033 


1.44 


1.03 


0.52 


0.52 


0.0002 


0.000 


0.425 


0.697 


0.923 


0.4 - 0.6 


2.47 


0.0090 


0.38 


1.57 


0.14 


0.27 


0.0122 


0.674 


0.675 


0.763 


0.851 


0.5 - 0.6 


2.16 


0.0083 


0.76 


1.31 


0.24 


0.28 


0.0084 


0.529 


0.615 


0.803 


0.935 


0.55 - 0.6 


1.87 


0.0094 


0.69 


1.27 


0.36 


0.36 


0.0002 


0.000 


0.575 


0.814 


0.979 


0.6 - 0.6 


1.54 


0.0032 


1.82 


1.03 


0.60 


0.60 


0.0002 


0.000 


0.461 


0.759 


1.006 


0.4 - 0.65 


2.87 


0.0088 


0.37 


1.67 


0.14 


0.28 


0.0121 


0.725 


0.728 


0.809 


0.895 


0.5 - 0.65 


2.50 


0.0082 


0.75 


1.38 


0.22 


0.32 


0.0090 


0.644 


0.682 


0.845 


0.971 


0.55 - 0.65 


2.29 


0.0087 


0.81 


1.32 


0.36 


0.36 


0.0012 


0.013 


0.630 


0.863 


1.017 


0.575 - 0.65 


2.06 


0.0092 


0.77 


1.28 


0.43 


0.43 


0.0002 


0.000 


0.604 


0.867 


1.039 


0.6 - 0.65 


2.02 


0.0082 


1.11 


1.20 


0.43 


0.43 


0.0002 


0.000 


0.607 


0.875 


1.060 


0.625 - 0.65 


1.99 


0.0089 


0.98 


1.22 


0.55 


0.55 


0.0002 


0.000 


0.594 


0.885 


1.082 


0.64 - 0.65 


1.94 


0.0014 


2.67 


1.04 


0.64 


0.64 


0.0002 


0.000 


0.535 


0.855 


1.096 


0.65 - 0.65 


1.82 


0.0031 


2.38 


1.04 


0.68 


0.68 


0.0002 


0.000 


0.503 


0.831 


1.095 


0.4 - 0.7 


3.22 


0.0087 


0.36 


1.75 


0.15 


0.29 


0.0118 


0.770 


0.780 


0.857 


0.943 


0.5 - 0.7 


2.96 


0.0078 


0.79 


1.45 


0.19 


0.35 


0.0091 


0.716 


0.752 


0.893 


1.012 


0.55 - 0.7 


2.76 


0.0077 


1.01 


1.35 


0.26 


0.36 


0.0081 


0.664 


0.724 


0.910 


1.051 


0.6 - 0.7 


2.51 


0.0078 


1.16 


1.27 


0.37 


0.38 


0.0035 


0.139 


0.687 


0.931 


1.097 


0.65 - 0.7 


2.28 


0.0084 


1.17 


1.23 


0.51 


0.51 


0.0002 


0.000 


0.653 


0.943 


1.144 


0.7 - 0.7 


2.04 


0.0030 


2.99 


1.03 


0.78 


0.78 


0.0002 


0.000 


0.542 


0.896 


1.177 


0.4 - 0.8 


4.05 


0.0082 


0.35 


1.93 


0.16 


0.31 


0.0116 


0.865 


0.888 


0.958 


1.045 


0.5 - 0.8 


4.18 


0.0074 


0.72 


1.68 


0.20 


0.40 


0.0103 


0.888 


0.896 


0.993 


1.102 


0.55 - 0.8 


4.27 


0.0073 


0.83 


1.64 


0.22 


0.43 


0.0098 


0.895 


0.904 


1.016 


1.136 


0.6 - 0.8 


4.09 


0.0069 


1.19 


1.50 


0.24 


0.45 


0.0081 


0.824 


0.882 


1.040 


1.182 


0.65 - 0.8 


3.72 


0.0066 


1.68 


1.36 


0.32 


0.47 


0.0071 


0.759 


0.844 


1.052 


1.214 


0.7 - 0.8 


3.34 


0.0074 


1.60 


1.32 


0.51 


0.51 


0.0028 


0.121 


0.776 


1.063 


1.257 


0.8-0.8 


2.49 


0.0029 


4.55 


1.03 


0.97 


0.97 


0.0002 


0.000 


0.601 


1.004 


1.334 


0.4 - 0.9 


5.36 


0.0076 


0.35 


2.15 


0.18 


0.33 


0.0112 


0.960 


0.987 


1.054 


1.142 


0.5 - 0.9 


5.17 


0.0070 


0.70 


1.83 


0.22 


0.43 


0.0099 


0.979 


1.007 


1.096 


1.204 


0.55 - 0.9 


5.30 


0.0067 


0.94 


1.74 


0.24 


0.47 


0.0093 


0.988 


1.013 


1.117 


1.236 


0.6 - 0.9 


5.49 


0.0065 


1.17 


1.69 


0.26 


0.53 


0.0088 


0.997 


1.011 


1.124 


1.255 


0.65 - 0.9 


5.71 


0.0061 


1.57 


1.61 


0.29 


0.55 


0.0083 


0.997 


1.020 


1.154 


1.297 


0.7 - 0.9 


5.15 


0.0059 


2.23 


1.45 


0.34 


0.58 


0.0068 


0.905 


0.981 


1.172 


1.337 


0.8 - 0.9 


4.48 


0.0065 


2.45 


1.35 


0.57 


0.61 


0.0042 


0.483 


0.907 


1.208 


1.417 


0.9 - 0.9 


3.19 


0.0012 


8.16 


1.01 


1.15 


1.15 


0.0002 


0.000 


0.715 


1.160 


1.517 


0.4 - 1.0 


6.81 


0.0069 


0.39 


2.34 


0.21 


0.37 


0.0107 


1.057 


1.081 


1.145 


1.233 


0.5 - 1.0 


6.55 


0.0066 


0.65 


2.06 


0.25 


0.46 


0.0094 


1.070 


1.106 


1.190 


1.302 


0.55 - 1.0 


6.60 


0.0063 


0.92 


1.91 


0.27 


0.52 


0.0086 


1.074 


1.115 


1.211 


1.333 


0.6-1.0 


6.59 


0.0060 


1.26 


1.79 


0.29 


0.59 


0.0082 


1.083 


1.121 


1.226 


1.356 


0.65 - 1.0 


7.11 


0.0056 


1.77 


1.72 


0.33 


0.65 


0.0077 


1.095 


1.116 


1.234 


1.377 


0.7- 1.0 


7.36 


0.0055 


2.03 


1.70 


0.36 


0.70 


0.0073 


1.099 


1.124 


1.254 


1.404 


0.8- 1.0 


7.06 


0.0050 


3.76 


1.48 


0.47 


0.76 


0.0056 


0.993 


1.095 


1.304 


1.484 


0.9 - 1.0 


6.34 


0.0048 


5.42 


1.33 


0.67 


0.83 


0.0049 


0.906 


1.057 


1.351 


1.570 


1.0-1.0 


4.97 


0.0065 


3.25 


1.34 


1.44 


1.44 


0.0001 


0.000 


0.928 


1.383 


1.674 


Varying Composition 


0.4 - 0.4 He 


0.46 


0.0043 


0.57 


2.36 


0.33 


0.33 


0.0003 


0.000 


0.307 


0.507 


0.671 


0.4 - 0.8 He 


2.92 


0.0090 


0.22 


3.89 


0.14 


0.29 


0.0123 


0.859 


0.890 


0.961 


1.050 


0.4 - 0.4 Mg 


1.15 


0.0026 


0.73 


0.74 


0.35 


0.35 


0.0003 


0.000 


0.319 


0.515 


0.674 


0.4 - 0.8 Mg 


4.47 


0.0080 


0.39 


1.54 


0.18 


0.33 


0.0112 


0.868 


0.885 


0.954 


1.040 


Varying oq 


0.6 -0.6 0.9a 


1.89 


0.0026 


2.01 


1.07 


0.60 


0.60 


0.0002 


0.000 


0.473 


0.775 


1.017 


0.6-0.6 1.1a 


1.23 


0.0002 


2.13 


0.95 


0.53 


0.53 


0.0002 


0.000 


0.509 


0.804 


1.023 


0.6 - 0.8 0.9a 


4.12 


0.0067 


1.42 


1.45 


0.27 


0.45 


0.0077 


0.802 


0.860 


1.037 


1.180 


0.6-0.8 1.1a 


3.98 


0.0071 


1.01 


1.54 


0.26 


0.44 


0.0093 


0.896 


0.906 


1.038 


1.173 


Synchronization 


0.6 - 0.6 synch 


1.32 


0.0104 


0.65 


1.16 


0.42 


0.42 


0.0002 


0.000 


0.552 


0.809 


1.001 


0.6 - 0.8 synch 


3.50 


0.0074 


0.95 


1.48 


0.31 


0.45 


0.0093 


0.896 


0.887 


1.027 


1.162 










Varying 


Viscosity 










0.6 - 0.8 lowvisc 


4.22 


0.0068 


1.20 


1.52 


0.18 


0.42 


0.0091 


0.881 


0.924 


1.063 


1.203 


0.6 - 0.8 highvisc 


4.45 


0.0069 


1.18 


1.55 


0.29 


0.44 


0.0084 


0.842 


0.877 


1.042 


1.186 










Varyin] 


g Numerical Resolution 










0.625 - 0.65 64k 


2.14 


0.0084 


1.11 


1.22 


0.49 


0.49 


0.0002 


0.000 


0.597 


0.879 


1.081 


0.625 - 0.65 128k 


1.98 


0.0084 


1.15 


1.18 


0.53 


0.53 


0.0002 


0.000 


0.574 


0.869 


1.084 


0.625 - 0.65 255k(2) 


1.92 


0.0003 


3.59 


0.99 


0.53 


0.53 


0.0002 


0.000 


0.587 


0.881 


1.084 


0.625 -0.65 510k 


2.10 


0.0003 


3.05 


1.04 


0.60 


0.60 


0.0002 


0.000 


0.556 


0.864 


1.085 
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Table? 

Table of merger remnant properties, continued. 



Sim 


(/)51 


(^disk)51 


(/ce)51 


(Alisk)49 


{^ce)49 






(P(f<'md))6 


(r(ra„d))8 


S^(ramd) 




0.4 - 0.4 


3.31 


2.41 


0.90 


6.13 


14.64 


0.0219 


0.0344 


0.00 


0.44 


0.08 


0.0071 


0.4 - 0.5 


6.40 


5.93 


0.47 


16.29 


8.83 


0.0168 


0.0292 


0.01 


0.48 


0.10 


0.0058 


0.5 - 0.5 


3.23 


2.29 


0.94 


7.37 


20.10 


0.0200 


0.0303 


0.00 


0.58 


0.11 


0.0059 


0.4 - 0.55 


7.47 


7.09 


0.38 


19.52 


7.69 


0.0154 


0.0302 


0.01 


0.46 


0.10 


0.0058 


0.5 - 0.55 


6.47 


5.84 


0.63 


14.77 


14.88 


0.0172 


0.0271 


0.01 


0.58 


0.12 


0.0051 


0.55 - 0.55 


3.26 


2.31 


0.95 


7.96 


23.00 


0.0191 


0.0280 


0.01 


0.67 


0.13 


0.0054 


0.4 - 0.6 


9.83 


9.57 


0.26 


24.33 


4.86 


0.0129 


0.0278 


0.01 


0.62 


0.11 


0.0060 


0.5 - 0.6 


6.52 


6.01 


0.51 


19.44 


12.68 


0.0159 


0.0278 


0.01 


0.66 


0.12 


0.0056 


0.55 - 0.6 


6.18 


5.57 


0.61 


17.01 


16.23 


0.0161 


0.0260 


0.02 


0.69 


0.14 


0.0050 


0.6 - 0.6 


3.50 


2.56 


0.94 


8.99 


25.58 


0.0181 


0.0260 


0.01 


0.75 


0.15 


0.0050 


0.4 - 0.65 


12.90 


12.66 


0.24 


26.67 


4.47 


0.0125 


0.0297 


0.01 


0.59 


0.10 


0.0062 


0.5 - 0.65 


7.42 


7.00 


0.42 


23.19 


11.15 


0.0145 


0.0282 


0.01 


0.72 


0.12 


0.0059 


0.55 - 0.65 


6.01 


5.52 


0.49 


21.47 


14.25 


0.0146 


0.0263 


0.02 


0.79 


0.14 


0.0053 


0.575 - 0.65 


5.92 


5.35 


0.57 


19.80 


16.65 


0.0153 


0.0257 


0.02 


0.82 


0.15 


0.0052 


0.6 - 0.65 


6.51 


5.87 


0.64 


17.81 


19.17 


0.0159 


0.0240 


0.02 


0.87 


0.16 


0.0047 


0.625 - 0.65 


5.05 


4.44 


0.61 


17.90 


19.32 


0.0151 


0.0236 


0.02 


0.87 


0.17 


0.0046 


0.64 - 0.65 


3.40 


2.62 


0.78 


13.31 


24.15 


0.0162 


0.0222 


0.03 


0.86 


0.18 


0.0042 


0.65 - 0.65 


2.90 


2.02 


0.88 


10.68 


27.20 


0.0168 


0.0237 


0.01 


0.92 


0.17 


0.0046 


0.4 - 0.7 


15.20 


14.97 


0.23 


28.60 


4.50 


0.0127 


0.0330 


0.01 


0.42 


0.09 


0.0057 


0.5 - 0.7 


8.41 


8.07 


0.34 


26.97 


9.52 


0.0133 


0.0263 


0.01 


0.82 


0.13 


0.0056 


0.55-0.7 


7.48 


7.05 


0.43 


24.88 


13.08 


0.0142 


0.0277 


0.01 


0.78 


0.13 


0.0057 


0.6 - 0.7 


6.15 


5.65 


0.50 


23.23 


16.20 


0.0144 


0.0259 


0.02 


0.85 


0.15 


0.0052 


0.65 - 0.7 


6.10 


5.51 


0.59 


20.09 


20.15 


0.0147 


0.0234 


0.03 


0.90 


0.17 


0.0045 


0.7 - 0.7 


2.77 


1.90 


0.87 


11.13 


30.06 


0.0161 


0.0225 


0.02 


1.06 


0.19 


0.0044 


0.4 - 0.8 


21.82 


21.63 


0.20 


33.03 


3.82 


0.0125 


0.0339 


0.01 


0.44 


0.09 


0.0058 


0.5 - 0.8 


12.33 


12.13 


0.20 


35.41 


5.32 


0.0106 


0.0277 


0.01 


0.76 


0.13 


0.0057 


0.55 - 0.8 


10.29 


10.09 


0.20 


36.41 


6.01 


0.0102 


0.0239 


0.02 


1.00 


0.16 


0.0052 


0.6 - 0.8 


7.77 


7.49 


0.29 


33.95 


10.18 


0.0115 


0.0233 


0.03 


1.04 


0.17 


0.0050 


0.65-0.8 


7.37 


6.98 


0.39 


30.15 


15.01 


0.0127 


0.0244 


0.02 


1.04 


0.17 


0.0050 


0.7 - 0.8 


6.88 


6.44 


0.44 


27.70 


18.42 


0.0126 


0.0238 


0.02 


1.08 


0.18 


0.0048 


0.8 - 0.8 


2.89 


2.04 


0.85 


11.66 


36.20 


0.0152 


0.0222 


0.01 


1.30 


0.21 


0.0044 


0.4 - 0.9 


30.16 


29.99 


0.17 


36.79 


3.55 


0.0126 


0.0356 


0.01 


0.42 


0.09 


0.0057 


0.5 - 0.9 


14.87 


14.70 


0.17 


40.15 


4.71 


0.0109 


0.0292 


0.02 


0.62 


0.12 


0.0052 


0.55 - 0.9 


10.82 


10.64 


0.17 


41.36 


5.44 


0.0101 


0.0262 


0.03 


0.69 


0.15 


0.0048 


0.6 - 0.9 


11.73 


11.56 


0.17 


42.56 


5.99 


0.0092 


0.0264 


0.02 


1.00 


0.15 


0.0056 


0.65 - 0.9 


9.50 


9.32 


0.18 


42.92 


6.92 


0.0089 


0.0213 


0.03 


1.33 


0.20 


0.0049 


0.7 - 0.9 


7.50 


7.24 


0.26 


38.80 


12.16 


0.0103 


0.0218 


0.03 


1.39 


0.20 


0.0049 


0.8 - 0.9 


6.88 


6.51 


0.37 


33.37 


19.44 


0.0111 


0.0208 


0.04 


1.40 


0.23 


0.0043 


0.9 - 0.9 


3.00 


2.31 


0.69 


16.81 


37.01 


0.0129 


0.0173 


0.08 


1.61 


0.31 


0.0033 


0.4 - 1.0 


38.03 


37.90 


0.13 


40.58 


3.04 


0.0109 


0.0385 


0.01 


0.40 


0.08 


0.0059 


0.5-1.0 


19.06 


18.92 


0.14 


44.81 


4.03 


0.0103 


0.0296 


0.02 


0.64 


0.13 


0.0052 


0.55-1.0 


13.40 


13.26 


0.14 


46.37 


4.58 


0.0096 


0.0274 


0.03 


0.75 


0.14 


0.0050 


0.6-1.0 


13.22 


13.08 


0.14 


47.85 


5.16 


0.0090 


0.0254 


0.04 


0.75 


0.16 


0.0045 


0.65- 1.0 


12.08 


11.93 


0.14 


48.14 


6.11 


0.0084 


0.0274 


0.03 


0.89 


0.15 


0.0050 


0.7-1.0 


11.13 


10.99 


0.14 


48.78 


6.75 


0.0080 


0.0237 


0.02 


1.52 


0.18 


0.0056 


0.8-1.0 


8.27 


8.06 


0.21 


44.97 


12.75 


0.0090 


0.0221 


0.03 


1.92 


0.21 


0.0054 


0.9-1.0 


7.03 


6.75 


0.27 


40.02 


18.84 


0.0094 


0.0190 


0.06 


1.96 


0.27 


0.0043 


1.0-1.0 


4.39 


4.04 


0.35 


31.75 


26.82 


0.0095 


0.0151 


0.17 


2.20 


0.39 


0.0033 


Varying Composition 


0.4 - 0.4 He 


3.89 


2.90 


1.00 


6.05 


15.22 


0.0230 


0.0367 


0.00 


0.30 


0.07 


0.0076 


0.4 - 0.8 He 


21.06 


20.86 


0.21 


34.21 


3.57 


0.0131 


0.0339 


0.01 


0.27 


0.09 


0.0057 


0.4 - 0.4 Mg 


3.28 


2.41 


0.86 


5.81 


14.66 


0.0215 


0.0339 


0.00 


0.47 


0.08 


0.0070 


0.4 - 0.8 Mg 


23.55 


23.36 


0.19 


32.37 


3.95 


0.0121 


0.0341 


0.01 


0.48 


0.09 


0.0057 












Varying 


ao 












0.6 - 0.6 0.9a 


2.38 


1.49 


0.89 


8.33 


24.59 


0.0177 


0.0235 


0.01 


0.87 


0.16 


0.0046 


0.6-0.6 1.1a 


3.93 


3.09 


0.84 


13.72 


22.58 


0.0176 


0.0245 


0.02 


0.64 


0.16 


0.0042 


0.6 - 0.8 0.9a 


9.01 


8.69 


0.33 


29.79 


11.99 


0.0121 


0.0238 


0.03 


0.94 


0.17 


0.0047 


0.6-0.8 1.1a 


10.76 


10.54 


0.22 


39.06 


7.16 


0.0102 


0.0243 


0.03 


0.83 


0.16 


0.0048 


Synchronization 


0.6 - 0.6 synch 


8.61 


7.94 


0.67 


18.95 


19.91 


0.0158 


0.0267 


0.01 


0.75 


0.14 


0.0052 


0.6 - 0.8 synch 


14.52 


14.30 


0.21 


40.21 


7.67 


0.0099 


0.0240 


0.03 


0.82 


0.17 


0.0047 


Varying Viscosity 


0.6 - 0.8 lowvisc 


5.03 


4.81 


0.22 


37.10 


7.18 


0.0102 


0.0206 


0.05 


1.19 


0.19 


0.0049 


0.6 - 0.8 highvisc 


6.89 


6.60 


0.29 


33.72 


10.43 


0.0115 


0.0232 


0.03 


1.01 


0.17 


0.0049 










Varyin 


I Numerical Resolution 










0.625 - 0.65 64k 


4.41 


3.79 


0.62 


17.59 


19.51 


0.0151 


0.0229 


0.02 


0.99 


0.17 


0.0050 


0.625 - 0.65 128k 


4.22 


3.56 


0.67 


16.24 


20.78 


0.0154 


0.0231 


0.02 


0.91 


0.17 


0.0047 


0.625 - 0.65 255k(2) 


4.77 


4.13 


0.63 


17.50 


19.73 


0.0153 


0.0234 


0.03 


0.84 


0.17 


0.0045 


0.625 -0.65 510k 


5.55 


4.85 


0.70 


15.92 


21.52 


0.0159 


0.0241 


0.02 


0.74 


0.17 


0.0043 
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